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ABSTRACT 

We extend color-magnitude relations for moderate-luminosity X-ray AGN hosts and non-AGN galaxies through 
the galaxy formation epoch iz ~ 1-4) in the Chandra Deep Field-North and Chandra Deep Field-South 
(CDF-N and CDF-S, respectively; jointly CDFs) surveys. This study was enabled by the deepest available 
X-ray data from the 2 Ms CDF surveys as well as complementary ultradeep multiwavelength data in these 
regions. We utilized analyses of color-magnitude diagrams (CMDs) to assess the role of moderate-luminosity 
ACjNs in galaxy evolution. First, we confirm some previous results and extend them to higher redshifts, finding, 
for example, that (1) there is no apparent color bimodality (i.e., the lack of an obvious red sequence and blue 
cloud) for AGN hosts from z w 0-2, but non-AGN galaxy color bimodality exists up to z « 3 and the relative 
fraction of red-sequence galaxies generally increases as the redshift decreases (consistent with a blue-to-red 
migration of galaxies); (2) most AGNs reside in massive hosts and the AGN fraction rises strongly toward 
higher stellar mass, up to z w 2-3; and (3) the colors of both AGN hosts and non-AGN galaxies become redder 
as the stellar mass increases, up to z w 2-3. Second, we point out that, in order to obtain a complete and reli- 
able picture, it is critical to use mass-matched samples to examine color-magnitude relations of AGN hosts and 
non-AGN galaxies. We show that for mass-matched samples up to z ~ 2-3, AGN hosts lie in the same region 
of the CMD as non-AGN galaxies; i.e., there is no specific clustering of AGN hosts in the CMD around the 
red sequence, the top of the blue cloud, or the green valley in between. The AGN fraction (« 10%) is mostly 
independent of host-galaxy color, providing an indication of the duty cycle of supermassive black hole growth 
in typical massive galaxies. These results are in contrast to those obtained with non-mass-matched samples 
where there is apparent AGN clustering in the CMD and the AGN fraction generally increases as the color 
becomes redder. We also find, for mass-matched samples, that the star-formation rates of AGN hosts are typi- 
cally a factor of « 2-3 larger than those of non-AGN galaxies at z ~ 0-1, whereas this difference diminishes at 
z ~ 1-3. With mass-selection effects taken into account, we find that almost all the results obtained in this work 
can be reasonably explained by two main ingredients, color-mass correlation (i.e.. X-ray AGNs preferentially 
reside in massive galaxies that generally tend to be redder than less-massive galaxies) and passive or secular 
evolution of galaxies. Our results show that the presence of moderate-luminosity AGN activity does not have 
a significant effect on the colors of galaxies and thus tightly constrain any effects from moderate-luminosity 
AGN feedback upon color-magnitude properties over the w 80% of cosmic time during which most of galaxy 
formation occurred. 

Subject headings: cosmology: observations — surveys — galaxies: normal — galaxies: active — galaxies: 
evolution — X-rays: galaxies 
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1. INTRODUCTION 

Deep multiwavelength surveys have greatly refined our un- 
derstanding of the cosmic star-formation and mass-assembly 
history. Measurements of galaxy luminosity functions at 
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multiple wavelengths have provided a reasonably tight con- 
straint on the evolution of the space density of the galaxy 
star-formation rate (SFR): the cosmic star formation history 
is known to within sa 30-50% up to redshifts of z ~ 1 and 
within a factor of « 3 at z « 1-6 (e.g., Hopkins 2004). It has 
been established that the galaxy SFR density has a broad max- 
imum at z ~ 1^, followed by a sharp drop from z ~ 1 to z ~ 
(e.g., Dickinson et al. 2003a; Giavalisco et al. 2004) and a 
fast decrease beyond z ~ 3^ (see, e.g., Hopkins & Beacom 
2006; Wilkins et al. 2008). It now seems clear that the bulk 
of present stellar mass was formed over a critical epoch at 
z « 1^ (e.g., Hopkins & Beacom 2006). 

Star-formation and active galactic nucleus (i.e., AGN) ac- 
tivity in galaxies are often found to go hand-in-hand (e.g.. 
Page et al. 2001; Alexander et al. 2005b; Netzer et al. 2007; 
Silverman et al. 2009). Both AGN activity and star forma- 
tion peak in about the same redshift range (z ~ 1.5- 2), and 
they undergo a similar decline below z ~ 1 (e.g., Cowie et 
al. 2003; Merloni et al. 2004; Silverman et al. 2008a). 
Mounting evidence has shown that supermassive black holes 
(SMBHs) centered in galaxies must play an important role 
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in galaxy evolution. The close connection between AGNs 
and their hosts is most strikingly shown by the tight corre- 
lations between the masses of SMBHs and the properties of 
their host-galaxy bulges (e.g., Haring & Rix 2004; Ferrarese 
& Ford 2005; Gultekin et al. 2009). Recent theoretical mod- 
els propose that AGN feedback is an important ingredient in 
understanding the interconnection between AGNs and their 
hosts (e.g., Silk & Rees 1998; Di Matteo et al. 2005; Bower 
et al. 2006; Croton et al. 2006; Fabian, Vasudevan, & Gandhi 
2008). Two main AGN feedback modes are thought to be 
significant. In the "radiative mode", the initially shrouded lu- 
minous AGN expels the obscuring gas (that feeds both the 
SMBH and star formation) and thus quenches star formation 
(e.g., Hopkins et al. 2005; Springel et al. 2005). The "ki- 
netic mode" invokes AGN heating to prevent hot gas from 
cooling and falling into a galaxy to form stars (e.g., Croton 
et al. 2006), which is generally thought to be important in 
the more massive halos and occurs at much smaller accretion 
rates than that of the radiative-mode feedback. To investigate 
the close connections between star formation and AGN activ- 
ity and thus fully understand galaxy formation and evolution, 
a sensible approach is to study both AGNs and their hosts. 

Galaxy colors are directly related to the star-formation, 
dust, and metal-enrichment history of galaxies and thus pro- 
vide important constraints for models of galaxy formation and 
evolution. The color-magnitude diagram (CMD; plot of rest- 
frame U -V/U -B colors vs. V-band/B-band absolute mag- 
nitudes) has proven to be one effective tool for exploring the 
role of AGNs in galaxy evolution. Previous CMD studies have 
obtained a number of results, e.g.: (1) The galaxy color bi- 
modality (i.e., separation of galaxies into the red sequence and 
the blue cloud) has been clearly seen in the CMD, while no 
color bimodality of X-ray AGN hosts appears to exist (e.g., 
Bohm & Wisotzki 2007; Nandra et al. 2007; Silverman et 
al. 2008b). (2) The clustering of AGNs in the CMD is dis- 
tinct. It has been shown that X-ray/optically-selected AGNs 
reside in the red sequence, the top of the blue cloud, and the 
green valley in between (e.g., Martin et al. 2007; Nandra et 
al. 2007; Rovilos & Georgantopoulos 2007; Westoby et al. 
2007); it has also been found that radio AGNs preferentially 
lie on the red sequence. X-ray AGNs lie in the green valley, 
and IR-selected AGNs are found in somewhat bluer hosts than 
X-ray AGNs (Hickox et al. 2009). (3) The broad distribution 
of host-galaxy colors of X-ray selected moderate-luminosity 
AGNs is dependent on the strong color evolution of lumi- 
nous (My < -20.7) galaxies and the influence of enhanced 
AGN activity in « 10 Mpc large-scale structures (Silverman 
et al. 2008b). (4) Morphologies of X-ray selected moderate- 
luminosity AGNs reveal no close connection between major 
mergers and AGN activity, and are consistent with minor in- 
teractions and/or secular evolution (e.g., Grogin et al. 2005; 
Georgakakis et al. 2008; Silverman et al. 2008b). (5) X-ray 
AGNs preferentially reside in luminous bulges (e.g., Nandra 
et al. 2007; Silverman et al. 2008b). However, we note 
that these CMD results are limited to redshifts below w 1 .4 
and relatively bright sources (usually mR < 24), primarily due 
to data-depth limitations. A natural follow-up question is 
what further insights can be gained by extending CMD stud- 
ies through the galaxy formation epoch when star-formation 
and AGN activity peaked? 

As data from exceptionally deep multiwavelength surveys 
have become available in the 2 Ms Chandra Deep Fields (see 
Brandt & Hasinger 2005 and Brandt & Alexander 2010 for 
reviews of deep extragalactic X-ray surveys), it is now feasi- 



ble to push CMD studies to higher redshifts. In this paper, we 
aim to assess the role of AGNs in galaxy evolution by means 
of CMD analyses through the galaxy formation epoch, using 
the 2 Ms CDF surveys (e.g., Alexander et al. 2003; Luo et al. 
2008) as well as superb complementary ultradeep multiwave- 
length data. A key aspect of this paper is the use of reliable 
stellar masses and the assessment of possible stellar-mass se- 
lection effects, because stellar mass is likely the most funda- 
mental observable parameter for understanding the properties 
of galaxies and there are suggestions of stellar-mass biases 
in CMD works (e.g., Silverman et al. 2009). This paper is 
structured as follows: § 2 describes the multiwavelength data 
used in this work; § 3 shows the derivations of source physi- 
cal properties as well as AGN identification; § 4 details sam- 
ple construction; § 5 gives the results obtained by this work, 
where we critically examine issues such as color bimodal- 
ity, AGN clustering in the CMD, and evolutionary trends in 
the colors of AGN hosts and non-AGN galaxies, with mass- 
selection effects taken into account; and finally, § 6 presents 
the conclusions and summary. Throughout this paper, all ab- 
solute magnitudes quoted are based upon the Vega magnitude 
system; X-ray luminosities are absorption-corrected (i.e., in- 
trinsic) and quoted in the 0.5-8 keV full band; and a cosmol- 
ogy of //o = 70.5 km s"' Mpc"', f^M = 0.274, and f^A = 0.726 
derived from the five-year WMAP observations (Komatsu et 
al. 2009) is adopted. 

2. MULTIWAVELENGTH DATA 
2.L Source Catalogs 

X-ray catalogs: We made use of the 2 Ms CDF point-source 
catalogs. The CDF-N catalog (Alexander et al. 2003) consists 
of 582 X-ray point sources (503 in the main Chandra source 
catalog and 79 in the supplementary optically bright Chan- 
dra source catalog). The CDF-S catalog (Luo et al. 2008) 
consists of 578 X-ray point sources (462 in the main Chandra 
source catalog, 86 in the supplementary CDF-S plus E-CDF-S 
Chandra source catalog, and 30 in the supplementary opti- 
cally bright Chandra source catalog). 

Optical/UV/IR catalogs: For the CDF-N, we used the 
Hawaii HDF-N optical and NIR catalog (U,B,V,R,I,z',HK'; 
Capak et al. 2004) as the base catalog that has a total of 
48,858 sources. This base catalog is sufficiently deep and 
complete to mR < 26. We cross-matched four other available 
optical/UV/IR catalogs to the Hawaii HDF-N catalog and kept 
only sources that have counterparts in the Hawaii HDF-N cat- 
alog. Those catalogs are (1) the GOODS-N ACS and IRAC 
catalogs (ACS F435W, F606W, F775W, F850LP, and IRAC 
3.6 ^m,4.5 fim, 5.8 fim, 8.0 ^m; Dickinson et al. 2003b);"''2 
(2) the GALEX HDF-N deep imaging survey catalog (NUV, 
FUV; GALEX Release 4 Data'^); (3) the ACS GOODS-N 
region /Q (< 24.5) catalog (Barger, Cowie, & Wang 2008); 
and (4) the GOODS-N MIPS 24 ^m catalog (Dickinson et 
al. 2003b). In this work we adopted the aperture-corrected 
photometry as detailed in the original references. 

For the CDF-S, we combined three catalogs to produce a 
base catalog: (1) the MUSYC By/?-detected optical catalog 
(U,B,V,R,I,o3,z; Gawiser et al. 2006); (2) the COMBO-17 
optical catalog (U,B,V,RJ + 12 medium-band filters; Wolf 
et al. 2004, 2008); and (3) the GOODS-S MUSIC catalog 
(U, F435W, F606W, F775W, F850LP, J,H,K„ and IRAC; 

" See http://archive.stsci.edu/pub/hlsp/goods/catalog_r2/. 
'^ See http://ssc.spitzer.caltech.edu/legacy/goodshistory.html. 
'3 See http://galex.stsci.edu/GR4/. 
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Grazian et al. 2006). All MUSYC sources and all the unique 
COMBO-17 and MUSIC sources (i.e., COMBO-17 and MU- 
SIC sources that did not match to a MUSYC source, with a 
matching radius of 0.5") were kept. With this approach, the 
base catalog has a total of 100,318 sources. This base cat- 
alog is also sufficiently deep and complete to wr < 26. We 
then cross-matched four other available catalogs to the base 
catalog and kept only sources that have counterparts in the 
base catalog. Those catalogs are (1) the MUSYC NIR catalog 
(J,H,K,; Taylor et al. 2009b); (2) the SIMPLE IRAC cata- 
log (Damen et al. 2010); (3) the GALEX CDF-S deep imag- 
ing survey catalog (NUV, FUV; GALEX Release 4 Data); 
and (4) the GOODS-S MIPS 24 ^m catalog (Dickinson et al. 
2003b). Again, we adopted the aperture-corrected photome- 
try for analysis. 

2.2. Cross-Matching 

Two cross-matching methods were used: closest- 
counterpart matching and likelihood-ratio matching. 
Closest-counterpart matching, as the name suggests, 
simply assigns the closest angular match as the counterpart 
(given some maximum matching radius); the separation 
between two sources is the sole criterion for this method 
(so it is easily implemented and executes quickly). Closest- 
counterpart matching performs acceptably (i.e., with a low 
false-match probability'"*), provided that an appropriate 
maximum-matching radius is adopted (see, e.g.. Table I). 
Here, closest-counterpart matching was used for the cross- 
matching between the base catalogs and other optical/UV/IR 
catalogs. 

The matching between the X-ray and optical/near-infrared 
sources is more challenging due to the fact that X-ray posi- 
tions are usually not as good as the optical/near-infrared ones, 
which are accurate to w 0.1-0.2"'^; the false-match probabil- 
ity increases because of the high density of faint background 
optical/near-infrared sources. Therefore, we used likelihood- 
ratio matching (e.g., Ciliegi et al. 2003; Luo et al. 2010) 
between the X-ray and base catalogs, which not only takes 
into account the positional accuracy, but also the expected 
magnitude distribution of the counterparts. Likelihood-ratio 
matching outperforms closest-counterpart matching for the 
matching between the X-ray and optical sources, especially 
when matching faint sources (see Table 1).'^ We note that the 
same matching approach was used for both the CDF-N and 
CDF-S, with similar results (shown in Table 1) produced for 

''^ We estimated the false-match probabihty by shifting the coordinates of 
base sources in RA and DEC by ±5" and re-correlating with sources from 
other catalogs. 

'^ As inferred from Table 1, X-ray positional accuracies are generally not 
as good as optical/IR ones, but are generally better than UV positions. We 
note that these X-ray positional accuracies are the best so far at these fluxes 
due to the long exposures of the 2 Ms CDFs: the majority of sources within 
the central /'encirded =8' area have a positional accuracy of < 0.4-0.5"; a 
number of sources even have a positional accuracy of < 0.2" that is compa- 
rable to the optical positional accuracy (e.g., Alexander et al. 2003; Luo et 
al. 2008). 

'* For the central ren^ircled =8' area of the CDF-N (CDF-S), we obtained 
a false-match probability of 13.78% (15.86%) and a median offset of 0.364" 
(0.291") using closest-counterpart matching (compared to 7.57% (7.04%) 
and 0.369" (0.302") using likelihood-ratio matching; see Table 1), with a 

matching radius of Tm = 1 -^y '^x '*' ^o '^^'^d, where Ax (Aq = 0. 1") is the 
X-ray (optical) positional error, and the coefficient of 1 .9 was chosen to obtain 
roughly the same number of matches as in the case of using likelihood-ratio 
matching for the purpose of direct comparison between methods. We refer 
readers to Luo et al. (2010) for more details on the likelihood-ratio matching 
method we used. 



Table 1 
Cross-Matching Results for the Central renckded ■ 
Areas 





Matching 


Median 


False 


Matching 


Catalog 


Radius ("; 


) Offset (") 


Rate 


Method 


(1) 


(2) 


(3) 


(4) 


(5) 


CDF-N 


GOODS-N ACS 


0.50 


0.145 


4.76% 


Closest-counterpart 


GOODS-N IRAC 


0.75 


0.202 


3.99% 


Closest-counterpart 


GALEX 


1.00 


0.505 


8.25% 


Closest-counterpart 


GOODS-N K, 


0.50 


0.202 


1.42% 


Closest-counterpart 


GOODS-N MIPS 


0.75 


0.276 


6.05% 


Closest-counterpart 


2 Ms CDF-N 


5.00 


0.369 


7.57% 


Likelihood-ratio 


CDF-S 


MUSYC NIR 


0.75 


0.167 


3.53% 


Closest-counterpart 


SIMPLE IRAC 


0.75 


0.215 


4.53% 


Closest-counterpart 


GALEX 


1.00 


0.475 


8.95% 


Closest-counterpart 


GOODS-S MIPS 


0.75 


0.338 


5.65% 


Closest-counterpart 


2 Ms CDF-S 


5.00 


0.302 


7.04% 


Likelihood-ratio 



All numbers quoted in this table are for matches to the base catalogs for 
sources with hir < 26 (only the central r^nckcled =8' radius areas around the 
respective average aim points of the 2 Ms CDFs are considered; see § 4.1 for 
source-selection areas). Columns: (1) Catalog that was matched to a base 
catalog. For the 2 Ms CDF catalogs, all X-ray sources (i.e., not limited to 
AGNs) were used for the assessments here. (2) Maximum matching radius 
(for closest-counterpart matching) or searching radius (for likelihood-ratio 
matching). (3) Median separation of all matches between a catalog and a 
base catalog. (4) False-matching probability. For the case of the 2 Ms CDF- 
N catalog, four filters (U.R,LHK') of the North base catalog were used with 
likelihood-ratio matching. The choice of these four filters was made to en- 
sure broad wavelength coverage and to cover all the sources in the North base 
catalog using as few filters as possible (this criterion of choosing filters also 
applies to the South base catalog). 7.57% is the average of the false rates 
derived with the four filters: 7.06% for J/-band, 8.66% for i?-band, 9.01% 
for /-band, and 5.57% for HK'-band, respectively. For the case of the 2 
Ms CDF-S catalog, eight filters (MUSYC B, MUSYC R, MUSIC /, MUSIC 
J, COMBO-17 U, COMBO-17 R, COMBO-17 915 nm, and IRAC 3.6/im) 
of the South base catalog were used with likelihood-ratio matching. 7.04% 
is the average of the false rates derived with the eight filters (in the afore- 
mentioned order, 8.61%, 9.28%, 6.83%, 5.88%, 7.31%, 7.04%, 7.31%, and 
4.05%, respectively). Note that the false rates from likelihood-ratio matching 
(i.e., 7.57% and 7.04%) not being smaller than those from closest-counterpart 
matching is due to the use of a larger maximum matching radius. (5) Method 
that was used for matching catalogs. 



each field. 



2.3. Redshift Catalogs 



We made use of both spectroscopic and photometric red- 
shifts. We preferred high-quality spectroscopic redshifts if 
available (see § 4.6.2 for more details); otherwise, we adopted 
photometric redshifts (Rafferty et al. 2010) derived with the 
Zurich Extragalactic Bayesian Redshift Analyzer (ZEBRA; 
Feldmann et al. 2006), using the ultradeep multiwavelength 
data described in § 2.1. ZEBRA has features to correct sys- 
tematic offsets in the photometry and construct new templates 
by modifying original templates based on the best fits between 
the photometry and original templates. We produced Zphot cat- 
alogs for both the CDF-N (using data from up to 18 filters) 
and CDF-S (using data from up to 42 filters). We note that 
our photometric redshifts have excellent quality considering 
dispersion and outliers (see Rafferty et al. 2010; also see 
§ 4.6.2). Luo et al. (2010) derived photometric redshifts for 
the 462 X-ray sources in the 2 Ms CDF-S main source cata- 
log in a similar way, but with a more sophisticated treatment 
of photometry; manual source deblending was performed and 
upper limits were set in all photometric bands without detec- 
tions. Therefore, photometric redshifts from Luo et al. (2010) 
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should on average have higher quality than our photometric 
redshifts and thus supersede our estimates if both estimates 
are available for a source (note that the two estimates are gen- 
erally in good agreement with each other). We refer readers 
to Rafferty et al. (2010) and Luo et al. (2010) for the details 
on Zphot derivation and a complete list of Zspec references used. 

3. SOURCE PROPERTIES 
3.1. Rest-Frame Absolute Magnitudes 

Using the data described in § 2.1 as well as upper limits in 
the K and IRAC bands when appropriate,'^ we computed rest- 
frame Johnson U,B,V absolute magnitudes (Mu,Mb,Mv) for 
each galaxy. We adopted the approach of template SED fit- 
ting to derive the absolute magnitudes, which was realized 
using ZEBRA. Compared to a linear (or log-linear) inter- 
polation/extrapolation method, the template SED fitting re- 
duces potential catastrophic failures, especially in cases of 
limited/incomplete spectral coverage. 

First, we constructed a comprehensive set of galaxy and 
galaxy/ AGN hybrid templates (see Rafferty et al. 2010 for 
details; also see Luo et al. 2010). Briefly, we used the 259 
PEGASE galaxy templates that were employed by Grazian 
et al. (2006), which cover different galaxy types (ellipti- 
cal, spiral, and starburst) and span a wide range of star- 
formation history and intrinsic extinction. We also used the 
10 AGN templates of Polletta et al. (2007) that include a 
variety of empirical quasar and Seyfert templates. To ac- 
count for both the nuclear and host-galaxy emission, we con- 
structed a series of hybrid templates between the represen- 
tative AGN templates and galaxy templates (e.g., Salvato et 
al. 2009). We chose 5 typical AGN templates (encompass- 
ing Type 1, Type 2, QSO 1, and QSO 2) and 16 typical PE- 
GASE galaxy templates (encompassing elliptical, spiral, and 
starburst) for the construction of hybrid templates. For each 
pair of AGN and galaxy templates, the templates were nor- 
malized by the total integrated fluxes, and four hybrids with 
varying galaxy/ AGN ratios (90:10, 75:25, 50:50, 25:75) were 
produced (note that the ratio of 10:90 is very similar to a 0: 100 
ratio, i.e., the original AGN templates, and thus was not in- 
cluded). Thus, we have a total of 330 galaxy/ AGN hybrid 
templates (i.e., 330=5 x 16 x 4+ 10); the number of final tem- 
plates is 259h-330=589. 

Second, we ran ZEBRA in Maximum-Likelihood mode 
to identify the best-fit template and compute absolute mag- 
nitudes for each source, with the input of photometry and 
redshift. For moderate-luminosity AGNs that are best fit by 
hybrid templates, we used only galaxy-component SEDs to 
estimate their host physical properties (e.g., absolute mag- 
nitudes). The difference between these property estimates 
and those derived using hybrid template SEDs (i.e., galaxy- 
component plus AGN-component) is either negligible or 
small because the optical and near-infrared emission from 
moderate-luminosity active galaxies is dominated by host 
starlight (see § 4.6.3 for details). Note that the absolute mag- 
nitudes output by ZEBRA are given in the AB system. In or- 
der to facilitate comparison with other works, these absolute 
AB magnitudes were converted into Vega magnitudes. 

'^ For each of the K and IRAC bands, a uniform upper hmit is applied 
if there is no detection at this specific band for a source. We find that the 
utilization of upper limits significantly helps constrain the best-fit template 
for sources with limited spectral coverage, and thus produces better estimates 
of physical properties (e.g., rest-frame absolute magnitudes, stellar masses, 
and SFRs) for these sources. 



We find our absolute magnitudes are in good agreement 
with those presented by Lehmer et al. (2008) for the over- 
lapping sources, who adopted the approach of convolving 
photometrically derived SEDs with Johnson U,B,V filter 
curves and directly computing rest-frame absolute magni- 
tudes through interpolations and/or extrapolations, rather than 
our approach of template SED fitting. When compared to 
Lehmer et al. (2008), for sources in the CDF-N (CDF-S), 
the median ratio between the two sets of absolute magnitudes 
is ss 1.0 (« 1.0) and there is a random scatter of < 0.2 mag 
(< 0.3 mag) [the median luminosity ratio is w 1.0 and ran- 
dom scatter is « 0.1 dex for both CDFs], which may be due 
to the fact that slightly different photometry and photometric 
redshifts were used. 

3.2. Stellar Masses 

We estimated a set of stellar masses at various bands (Mx.*) 
for each source using the following tight correlations between 
rest-frame optical colors and stellar mass-to-light ratios, 

log(MA.,/MQ) = log(LA/iA.0) + ^A(MB-Mv) + flA-O.lO, 

(1) 

where X = B, V, R, I, J, H, and /T-band, respectively. Equa- 
tion (1) was prescribed in Table Bl of Zibetti, Chariot, & Rix 
(2009) where the values of the coefficients (ax and bx) can 
be found. The normalization in Eqn. (1) has been adjusted by 
-0.10 dex to account for our adopted Kroupa (2001) initial 
mass function (IMF). We further used the relation Lx/Lx^q = 
2.512*^^ Q"*^^ for the mass estimation. We note that Table Bl 
of Zibetti, Chariot, & Rix (2009) is a direct update of Ta- 
ble 7 of Bell et al. (2003) and was derived using a technique 
that constructs spatially resolved maps of stellar-mass surface 
density in local galaxies based on optical and NIR imaging, 
where the latest stellar population synthesis models were in- 
corporated. We chose to use Table Bl of Zibetti, Chariot, & 
Rix (2009) rather than Table 7 of Bell et al. (2003) because 
we find, when using the prescription in Zibetti, Chariot, & Rix 
(2009), that (1) for the same color (e.g., B-V in Eqn. (1)), the 
stellar-mass estimates derived at different bands are in better 
agreement with each other, in terms of both scaling and tight- 
ness; and (2) for different colors (i.e., B-V and B-R; the 
set of the coefficients, ax and bx, for the B-R color can also 
be found in Table Bl of Zibetti, Chariot, & Rix 2009)l^ die 
stellar-mass estimates derived at the same band are in better 
agreement with each other, in terms of scaling. 

Studies have shown that /T-band galaxy luminosities are 
5-10 times less sensitive to dust and stellar-population effects 
than optical luminosities (e.g.. Bell & de Jong 2000), which 
allows accurate stellar-mass estimations for galaxies. There- 
fore, for each source, we consistently adopted the stellar mass 
derived with the reddest absolute magnitude as the final mass 
estimate, applying the condition that the source must actually 
have a corresponding rest-frame detection (i.e., Mk,* has the 
highest priority and is adopted if a source has a rest-frame 
TT-band detection). Upper limits in the K and IRAC bands 
used in the template SED fitting process (see § 3.1) were not 
counted for the purpose of determination of reddest rest-frame 
coverage, although they do help avoid unrealistically large 
mass estimates for the sources with limited spectral coverage. 

'* We also estimated another set of stellar masses at various bands for each 
source using the B-R color, again according to Table Bl of Zibetti, Chariot, 
& Rix (2009). We find that the agreement between the two sets of stellar- 
mass estimates is good, so we adopted the B — V color, i.e., Eqn. (1), for the 
stellar-mass estimates, given the fact that the rest-frame V-band coverage is 
better than the rest-frame R-band coverage for the sources. 
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For sources in the CDF-S, we compared our stellar-mass 
estimates with those presented by Borch et al. (2006) and find 
general agreement between methods, with a median ratio of 
« 0.84 between our and their estimates and < 0.4 dex random 

scatter, which is adequate for our purposes in this work.'^ We 
speculate that the offset of mass scale may arise from the sub- 
tle differences in the IMF and the random scatter may mostly 
be due to the fact that we are using different photometric red- 
shifts from those used by Borch et al. (2006). 

3.3. Star-Formation Rates 

We estimated the star-formation rate (SFR) of each source 
following Eqn. (1) in Bell et al. (2005): 

SFR(Mo yr-i) = 9.8 x IO-'HLuy/Lq+Lir/Lq), (2) 



where the solar bolometric luminosity L0 = 3 .9 x 1 0^^ 



ergs 



The dust-uncorrected ultraviolet luminosity Luv was com- 
puted following Luv = 3.3iylA2^0Q A) (see §3.2 of Bell et al. 

2005), where /^(2800 A) is the rest-frame 2800 A monochro- 
matic luminosity that was estimated using the best-fit tem- 
plate SEDs derived in § 3.1. For the estimation of the 
8-1000 /im infrared luminosity Lir, we employed an IDL 
routine "chary _elbaz_24um.pro"^" provided by Chary & El- 
baz (2001), which incorporates a library of 105 template 
SEDs that reproduce the observed 24 /im luminosity-LiR cor- 
relations for local galaxies. We utilized observed MIPS 
24 /im fluxes (or upper limits) and redshifts to derive Lir 
(or upper limits on Lir). Papovich et al. (2007) found that 
SFRs derived with Eqn. (2) (i.e., SFRuv+24/im) are overesti- 
mated for luminous IR galaxies, compared to the estimates of 
SFRuv+24,70,i60/im that also consider longer- wave length MIPS 
bands (i.e., 70 and 160 /im). They corrected the overestimates 
using an empirical second-order polynomial that fits the corre- 
lation between log(SFRuv+24,70,i60Mm) and log(SFRuv+24A<m)- 
We therefore applied the correction derived by Papovich et al. 
(2007) to our SFRuv+24/im estimates that are > 25OM0 yr"'. 
A few percent of the sources were corrected downward, with 
a multiplicative correction factor ranging from ^0.3-1.0 (the 
median correction factor is ~ 0.8). 

For sources with MIPS 24 /im detections in the CDF-N 
(CDF-S), we find good agreement between our SFR values 
and those derived by Lehmer et al. (2008), with a median ra- 
tio of « 1.0 (« 1.0) between our and their SFRs, and < 0.2 
(< 0.3) dex random scatter. For sources without MIPS 24 /im 
detections in the CDF-N (CDF-S), we find « 100% (« 100%) 
of our upper limits on SFR are consistent with those presented 
by Lehmer et al. (2008). As a further check, we also com- 
puted SFRuv+70/^m (i-c, using observed MIPS 70 /im fluxes 
to derive Lir; an IDL routine "chary _elbaz.pro"^" was used) 
for a subsample of 475 sources in the CDF-S base catalog 
that have both MIPS 24 /im and 70 /im detections, and com- 
pared with their corresponding SFRuv+24;^m values. We find 
good agreement between these two sets of SFR estimates; 
the median ratio between SFRuv+vo/im and SFRuv+24/jm is 
« 1.07, with a random scatter of < 0.4 dex and no system- 
atic deviations. This agreement further demonstrates that our 

" We note that estimates of stellar mass using photometry are generally 
accurate to within a factor of ss 3-5. We refer readers to other works for de- 
tailed discussion of the complexity and uncertainty of stellar-mass estimates 
(e.g., Bundy et al. 2006; Gallazzi & Bell 2009; Muzzin et al. 2009; Zibetti, 
Chariot, & Rix 2009). 

^'' See http://david.elbaz3.free.fr/astro_codes/chary_elbaz.html. 



SFRuv+24pm estimates are rehable and adequate for our pur- 
poses in this work. 

3.4. AGN Identification 

For the 2 Ms CDF point-source catalogs, we used four pri- 
mary criteria that rely upon distinct AGN physical properties 
to identify AGN candidates. The techniques of AGN iden- 
tification used here are detailed in, e.g., Bauer et al. (2004) 
and Lehmer et al. (2008). Here, we only describe these tech- 
niques briefly. We note that a source can be identified as an 
AGN through multiple criteria. 

(a) X-ray Luminosity: The intrinsic X-ray luminosity (Lx; 
quoted in the 0.5-8 keV band throughout this paper) was es- 
timated using the following equation: 



L0.5-& keV = 47r(i^/o.5-8 keV,int( 1 + z) , 



(3) 



where di is the luminosity distance, /0.5-8 kev.int is the intrinsic 
(i.e., absorption-corrected) flux, and F is the power-law pho- 
ton index of the X-ray spectrum. We derived /0.5-8 kev.int fol- 
lowing these procedures: (1) Using XSPEC (Arnaud 1996), 
we modeled the X-ray emission with an absorbed power-law 
model (both Galactic and intrinsic absorption), which is de- 
fined as zpowxwabsxzwabs in XSPEC. The photon index 
was fixed to F = 1.8 (typical for intrinsic AGN spectra) for 
each source and the redshifts of the zpow and zwabs compo- 
nents were fixed to that of the source. We additionally fixed 
the Galactic column density to A^h = 1.6 x 10^*^ cm"^ for the 
CDF-N and to A^h = 8.8 x 10''' cm'^ for the CDF-S (Stark et 
al. 1992). (2) We used the above absorbed power-law model 
to find the intrinsic column density that reproduces the ob- 
served band ratio, which is defined as the ratio of count rates 
between the hard (2-8 keV) and soft (0.5-2 keV) bands. (3) 
We removed both Galactic and intrinsic absorption to obtain 
the intrinsic flux. We find our /0.5-8 keV.int estimates to be in 
good agreement with those derived from spectral fitting utiliz- 
ing ACIS Extract^' (F. Bauer et al., in preparation), in terms 
of correlation and dispersion. 

Local X-ray observations show that purely star- forming 
non- AGN galaxies usually do not have intrinsic Lx > 3 x 10"*^ 
erg s"'. We thus classify a source with Lx > 3 x lO'*^ erg s"' 
as an AGN (i.e., a luminous AGN). We note that there could 
be potential intruders such as high-redshift strongly star- 
forming sources (e.g., submillimeter galaxies) that have a sub- 
stantial amount of X-ray emission mainly from active stel- 
lar populations (e.g., Alexander et al. 2005a). However, the 
Lx-to-SFR ratio of these star-forming galaxies is typically 
lower than that of AGNs (see Criterion (c) below), which 
can be used as a discriminator to differentiate these non- AGN 
galaxies from AGNs. Indeed, we find that over « 92% of 
the 2 Ms CDF AGNs identified with Criterion (a) were also 
identified with the Lx-to-SFR ratio criterion (i.e.. Criterion 
(c)), which suggests that contamination by high-redshift star- 
forming galaxies with Criterion (a) is not significant in our 
case. 

(h) X-ray Hardness: One signature of moderately-to-highly 
obscured (A^h > 10^^ cm"^) AGNs is their hard X-ray spectra. 
If a source has an effective photon index of Feff < 1 .0, we flag 
it as an AGN (i.e., an obscured AGN). 

(c) X-ray-to-SFR Correlation: AGNs often show signifi- 
cant X-ray luminosity excesses over what is expected based 

^' See http://www.astro.psu.edu/xray/docs/TARA/ae_users_guide.html 
for details on ACIS Extract (Broos et al. 2010). 
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on the Lx-SFR correlation. We modified the Lx-SFR corre- 
lation shown in Persic & Rephaeli (2007) to Lo.5-8kev/SFR = 
1.06 X lO'**' erg s"' (Mq yr"')"', taking into account the fact 
that we computed SFRs using Eqn. (2). If a source has 
io.5-8kev/SFR > 3 times that from the above Lx-SFR corre- 
lation, which corresponds to a factor of > 2.5 times the RMS 
scatter of the Persic & Rephaeli Lx-SFR correlation, we clas- 
sify it as an AGN. 

(d) X-ray-to-Optical Flux Ratio: The X-ray-to-optical 
flux ratio log(/o.5-8 keV,int//R) = "1 IS regarded as a useful 
AGN/galaxy discriminator. We mark a source as an AGN if it 
satisfies log(/o.5-8keV,int//R) > -1- 

4. SAMPLE CONSTRUCTION 

We constructed a parent sample within the 2 Ms CDFs that 
consists of 207 X-ray selected moderate-luminosity AGNs 
and a population of 12,329 galaxies that the AGNs were 
drawn from. We then extracted four samples, including both 
stellar-mass and luminosity selected samples, from this parent 
sample for various purposes. The details of sample construc- 
tion are described below. 

4.1. Source-Selection Areas 

We restricted our source selection to the central rencircied =8' 
radius areas around the respective average aim points of the 
2 Ms CDFs that are weighted by exposure time (for the 
CDF-N: aj2ooo.o = 12h36'"45.^7, Jj2ooo.o = +62°13'58"; for 
the CDF-S: ajzooo.o = 03h32™28.^8, <5j2ooo.o = -27°48'23"), 
and we also required GOODS-N/GOODS-S coverage with 
HST of sources (only k, 15% of the central rencircied =8' ar- 
eas do not have GOODS coverage). The corresponding to- 
tal selected area is ss 314 arcmin^ (w 155 arcmin^ for the 
CDF-N and « 159 arcmin^ for the CDF-S). The reasons 
for this areal choice are (1) The central CDF areas provide 
the deepest X-ray coverage to date, with sensitivity limits of 
< 2.3 X 10"'^ erg cm"^ s"' for the 0.5-8 keV band. Such deep 
X-ray data can reduce the bias that shallower X-ray surveys 
have against finding AGNs in less-massive hosts (e.g., Mer- 
loni & Heinz 2008). X-ray surveys have been extremely effi- 
cient in defining reliable AGN samples; the 2 Ms CDFs have 
reached an AGN density of ss 7200 deg~^ (e-g-, Bauer et al. 
2004). (2) More accurate X-ray positions for sources in the 
central areas lead to more reliable matching between X-ray 
and faint optical sources (see § 2.2). (3) The redshift success 
rate (i.e., availability of spectroscopic redshifts, Zspec) is rela- 
tively high due to intensive spectroscopic investment in these 
areas. Furthermore, the quality of photometric redshifts (Zphot) 
has been greatly improved utilizing the superb multiwave- 
length data available. (4) The deep multiwavelength coverage 
within the GOODS areas (e.g., ultradeep HST observations 
provide the best optical photometry and images for analysis) 
allows for the construction of well-sampled rest-frame spec- 
tral energy distributions (SEDs) that are key to estimations of 
Zphot and physical properties. 

For the CDF-N, there are 11,951 galaxies (253 AGNs) 
located within the northern rgncircied =8' area and having 
GOODS coverage; for the CDF-S, there are 21,749 galaxies 
(250 AGNs) located within the southern rencircied =8' area and 
having GOODS coverage. ^^ 

^^ Note that the total galaxy number in the respective central r„,circled =8' 
area and having GOODS coverage significantly differs between the North 
and South base catalogs, i.e., 11,951 vs. 21,749, which is mainly due to the 



4.2. Filter-Number Cut 

To ensure the reliability of the observed source SEDs, we 
discarded sources that have detections in less than five fil- 
ters. Consequently, the above galaxy numbers of 11,951 
and 21,749 are reduced to 1 1,949 (« 0.02% loss) and 21,376 
(ss 1.72% loss), respectively. We do not expect this filter- 
number cut to introduce any significant bias because only very 
small fractions of sources were dropped, and dropped sources 
are typically very faint (and will be excluded later anyway; 
see the m^ < 26 cut in § 4.3). We note that the numbers of 
AGNs in the CDF-N and CDF-S are not affected by this cut 
(still 253 and 250, respectively). 

4.3. Redshift Cut and Magnitude Cut 

In addition to our areal choice and filter-number cut, we 
limited our sample selection to a redshift range of < z < 4. 
We further applied an /?-band magnitude cut of niR < 26 for 
the source selection, which ensures that the sources in our 
sample have high-quality photometry and photometric red- 
shifts, and are only subject to minimal cross-matching fail- 
ures. Complexity in the application of a uniform cut of 
mR < 26 arose from the fact that different /?-band filters were 
used for the CDF-N (the Subaru R broad-band filter) and 
CDF-S (the COMBO- 17 and MUSYC R broad-band filters) 
observations. We converted the two CDF-S /?-band magni- 
tudes to the Subaru /?-band magnitude using a /T-correction 
package (kcorrect.v4_l_4^^; Blanton & Roweis 2007), which 
convolves the photometrically derived SED (see §2.1 for the 
photometry catalogs used) of a source with the above three 
/?-band filter curves and computes the difference (Awr) be- 
tween the derived /?-band magnitudes. Typically, Am^ is 
small (|AmR| < 0.15) with a median value of « 0.02. After 
applying the redshift and magnitude cuts, our sample consists 
of 401 AGNs and 19,202 galaxies. 

4.4. Stellar-Mass Cut and/or Luminosity Cut 

Stellar mass is likely the most fundamental observable pa- 
rameter for understanding the properties of galaxies. Us- 
ing stellar-mass selected samples to study galaxy evolution 
has some advantages over using either color- or luminosity- 
selected samples, e.g., avoiding biases associated with sam- 
ple selections (see § 6 for details). Figure 1 shows stellar 
mass as a function of redshift for our sample of 401 AGNs and 
their 19,202 parent galaxies (see § 4.3). According to Fig. 1, 
our sample of blue (red) galaxies is roughly complete above 
lO'^-^M© (IO'^-^Mq) from z « to z « 4 (red and blue galax- 
ies are separated using Eqn. (4) in § 5.1.1). The completeness 
limit of lO'^'^M© (10 M©) is confirmed using the technique 
presented in § 5.1 of Meneux et al. (2008) that constructs the 
mass completeness limit as a function of redshift based on the 
properties (i.e., rest-frame absolute magnitude, redshift, and 
stellar mass) of galaxies in a sample. We therefore applied a 
stellar-mass cut of > IO^'^Mq to our sample. 

To make direct comparisons with previous CMD studies 
that typically made use of luminosity limited samples, we also 
kept sources that satisfy a luminosity cut of My < -19 (i.e., 
we included both stellar-mass and luminosity selected sam- 
ples in our parent sample; see § 4.6.1). Therefore, there are 

fact that the South base catalog was constructed in a different way from the 
North, and it includes a large number of faint sources (e.g., mir > 26) because 
the GOODS-S MUSIC catalog has a set of fainter optical and near-infrared 
detection limits than the North base catalog. 
^' See http://cosmo.nyu.edu/blanton/kcorrect/. 
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Fig. 1 . — Stellar mass as a function of redshift for the sample that consists 
of 401 AGNs and their 19,202 parent galaxies (see § 4.3). AGNs are shown 
as black large filled circles; red/blue galaxies are shown as red circles/blue 
dots (red and blue galaxies are separated using Eqn. (4) in § 5.1.1). The solid 
(dashed) fine indicates the stellar-mass cut of 10''^Mq (IO'^'^Mq), above 
which blue (red) galaxies in Sample A are roughly complete (see § 4.6). 
The apparent redshift quantization reflects the logarithmic redshift steps we 
adopted for deriving the photometric redshifts. (A color version of this figure 
is available in the onhne journal.) 



12,329 galaxies (375 AGNs) in our sample, which have either 
M* > IO'^-^Mq orMv < -19. 

4.5. X-ray Luminosity Cut and Removal of Broad-Line AGNs 

We further chose to restrict the luminosity range of AGNs 
to 41.9 < log(Lx/(erg s"')) < 43.7 (Silverman et al. 2008b). 
The choice of the lower and upper luminosity limits was 
made to remove the contamination of any non-AGN galax- 
ies having significant X-ray emission, to minimize potential 
luminosity-dependent effects, and to ensure that the optical 
and near-infrared emission is primarily due to host galaxies 
(see § 2.2 of Silverman et al. 2008b for details; also see 
§ 4.6.3). There are 215 AGNs (out of 375; see § 4.4) that 
satisfy this X-ray luminosity cut: 215 = 375-89-71, where 
89 AGNs have log(Lx/(ergs-')) < 41.9 and 71 AGNs have 
log(Lx/(erg s-i)) > 43.7. We note that all of the 89 AGNs 
with log(Lx/(erg s"')) < 41 .9 dropped have z < 1 .0, which is 
not the redshift range of primary interest for this study. Fur- 
thermore, we have verified that our basic results do not change 
with these 89 AGNs included. 

Even with all the above source-filtering procedures, a few 
broad-line AGNs (BLAGNs) can still survive in our sample 
that are optically bright and significantly contribute to the 
optical emission; we cannot obtain reliable stellar-mass con- 
straints for these sources. We thus used the identifications 
of BLAGNs in the CDFs (Bai-ger et al. 2003; Szololy et al. 
2004) to remove the BLAGNs remaining in our sample. After 
this screening, there are 207 AGNs remaining in our sample 
(i.e., 8 BLAGNs were dropped). 

4.6. Sample Properties 

4.6.1. Sample Definitions 

Using the above criteria (/"enciicied = 8', having GOODS cov- 
erage, having no less than 5 filters, < z < 4, wr < 26, 
M* > IO'^-^Mq orMv < -19, 41 .9 < log(Lx/(erg s"')) < 43.7, 
and not being a BLAGN; the latter two criteria are applied 
only to AGNs), we constructed a parent sample that consists 



of 207 X-ray selected AGNs, which was drawn from a popu- 
lation of 12,329 galaxies. Table 2 shows a summary of parent 
sample construction. 



Table 2 
Summary of Parent Sample Construction 





Galaxies 


Galaxies 


AGNs AGNs 


Criteria 


(CDF-N) (CDF-S) (CDF-N) (CDF-S) 


I. ''encircled = 8', GOODS Coverage 


11,951 


21,749 


253 250 


2. > 5 filters 


11,949 


21,376 


253 250 


3. < z < 4, mR < 26 


8588 


10,614 


230 171 


4. M*> lO'-5M0orMv < -19 


5738 


6591 


214 161 


5. 10^' « < Lx < 10« ' 


— 


— 


116 99 


6. Not a BLAGN 


— 


— 


115 92 


Parent Sample 


5738 


6591 


115 92 



We extracted four samples (A, B, C, and D) from the par- 
ent sample that are appropriate for various purposes. Table 3 
shows some basic information about these four samples. In 
Sample A (0 < z < 4 and M^ > lO^-^M©), blue galaxies with 
M-i, > 1O'^M0 are roughly complete from z ~ 0^, whereas 
the completeness of red galaxies is not as high (see Fig. 1). 
We used Sample A mainly for the examination of possible 
galaxy color bimodality up to high redshifts (see § 5.1) be- 
cause Sample A, consisting of a sufficient number of blue and 
red galaxies, is suitable for this purpose. There are only a 
few AGNs with 3 < z < 4 in Sample A. We thus focused on 
a redshift subrange of < z < 3 in Sample A for studying 
both AGNs and their hosts, using Sample B (0 < z < 3 and 
M* > IO'^-^Mq) and Sample C (0 < z < 3 andM,, > 10^°-^Mq) 
that are two subsets of Sample A. We note that the mass cut 
of IO^'-'Mq in Sample C is the mass above which both AGN 
hosts and non-AGN galaxies have similar distributions of rest- 
frame U-V color (see § 5.3). In Sample B, only blue galax- 
ies are roughly complete; while in Sample C, both blue and 
red galaxies are roughly complete (see Fig. 1). Therefore, 
these two samples were often used together to examine the 
effects of sample incompleteness on various results (see § 5.2 
and 5.3). Sample D (0 < z < 3 and My < -19) was used to 
represent luminosity limited samples that have typically been 
adopted in previous CMD studies; the CMD results obtained 
with Sample D were compared with those from mass-limited 
samples (i.e. Samples A, B, and C, which are the main fo- 
cus of this work). In the remainder of § 4.6, we describe the 
properties of Sample A, which includes both Sample B and 
Sample C. 

4.6.2. Quality of Photometric Redshifts 

We assessed the Zphot quality for Sample A (see Fig. 2) by 
using a number of quantities produced by a comparison of 
the photometric and spectroscopic redshifts of the sample: 
the normalized median absolute deviation (unmad = 1.48 x 
median[|(Az-median(Az))/(l+Zspec)|]), the average absolute 
scatter (AAS = mean[| Az/(1 +Zspec)|]), and the percentages of 
outUers (P().2(0.i)_outiiei-) with |Az|/(l+Zspec) > 0.2(0.1), where 
Az = Zphot-Zspec- For the AGNs, we find that ctnmad = 0.0154, 
AAS = 0.0214, Po.2_outiier = 1 -35%, and Po.i_outiier = 4.05%; for 
the galaxies, we find that (Tnmad = 0.0153, AAS = 0.0323, 
A).2_outiier = 2.54%, and Po.i_outiiei- = 5.51%. We owe such a 
good quality of Zphot to the large CDF multiwavelength obser- 
vational investments and continuous efforts of improving Zphot 
estimation techniques (see, e.g., Luo et al. 2010; Rafferty et 
al. 2010; and references therein). However, as cautioned by 
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Table 3 
Galaxy and AGN Samples Drawn from Parent Sample 



Sample 



Condition 



No. of Galaxies No. of AGNs 



Main purpose 



Sample A < z < 4, M* > lO'^-^M© 

Sample B < z < 3, M* > lO^-^M© 

Sample C 0<z<3,M^> W^'^-^Mq 
Sample D < z < 3, My < -19 



4357 


185 


1468 


139 


11,119 


202 



Color bimodality 
CMD, color dependence, AGN fraction, etc. 
CMD, color dependence, AGN fraction, etc. 
CMD, color dependence, AGN fraction, etc. 



Luo et al. (2010) and Rafferty et al. (2010), the spectroscopic 
subsample is likely to have significantly better photometric 
redshifts than the full sample due to two reasons: (1) the 
sources in the spectroscopic subsample are generally brighter 
and therefore not entirely representative of the full sample; 
and (2) the template-improvement step (see § 2.3) used in the 
Zphot derivation optimizes the templates for the spectroscopic 
subsample and thus introduces a bias that is favorable to the 
spectroscopic subsample. ^"^ 

We do not expect the photometric redshift uncertainties to 
introduce a systematic bias in the estimation of source proper- 
ties (e.g., rest-frame absolute magnitudes) because there is no 
systematic offset between photometric redshifts and spectro- 
scopic redshifts and the residuals between photometric red- 
shifts and spectroscopic redshifts are approximately symmet- 
ric (see Fig. 2). Table 4 shows Zspec fractions of the sources 
within different redshift ranges in Sample A, which indicates 
that the Zspec fractions for AGNs are quite high and those for 
galaxies are reasonably high, especially at low redshifts. 



Table 4 
Zspec Fractions for Different Redshift Bins in Sample A 



Redshift No. 


of Galaxies (% of 


Zspec) No. 


ofAGNs(% of Zspec) 


0<z< 1 


1626 (68.7%) 




83 (94.0%) 


1<Z<2 


1570 (36.5%) 




85 (71.8%) 


2<z<3 


1161 (15.0%) 




17 (47.1%) 


3<z<4 


584(10.6%) 




3 (33.3%) 


0<z<3 


4357 (42.8%) 




185 (79.5%) 


0<z<4 


4941 (39.0%) 




188 (78.7%) 



4.6.3. AGN Contamination 

The X-ray luminosity cut of 41.9 < log(Lx) < 43.7 (see 
§ 4.5) was adopted to minimize potential AGN contamina- 
tion of optical and near-infrared emission from host galax- 
ies. Below we present four additional arguments to show 
that the optical and near-infrared emission from the 188 ac- 
tive galaxies in Sample A is dominated by host starlight and 
not strongly influenced by the AGNs. (1) Of these 188 AGNs, 
186 (« 99%) have best-fit SED templates with the optical and 
near-infrared emission dominated by galaxy starlight. (2) Fig- 
ure 3 shows 0.5-8 keV luminosities of these 188 AGNs as 

^^ In order to assess the actual quality of our photometric redshifts, we 
repeated a divided-sample test several times: we ran template-improvement 
mode with fs 3/4 of the sources in the spectroscopic subsample that were 
randomly selected and then applied the obtained improved templates to the 
unselected ks 1 /4 of the sources to test for quality. The tests suggest that the 
percentage of outliers for the non-trained sources is larger than that for the 
trained sources: Pq 2_outlier ~ 6.5% (ss 14.9%) for the non-trained galaxies 
(AGNs). 



a function of redshift, V-band absolute magnitude, and rest- 
frame U -V color (i.e., Mu-My), respectively. According to 
Figs. 3(b, c), there is no clear correlation between X-ray lu- 
minosity and either My or L'^-V. This is consistent with the 
X-ray and optical emission not being produced by the same 
process, suggesting that the optical emission is not likely to 
be dominated by the AGN component. (3) Following Silver- 
man et al. (2008), we determined conservative upper limits on 
the AGN contribution to the total (galaxyn-AGN) optical emis- 
sion using the HST ACS Veoe- and zgso-band images (rescaled 
to 0.03" pixel"') provided by GOODS. About a dozen AGNs 
were not included in the analysis due to their proximity to 
the ACS field edge. We measured optical counts in circu- 
lar apertures of two different sizes (r = 3 and 25 pixels, i.e., 
r = 0.09" and 0.75", respectively) located at the centroid of 
the optical emission, r = 0.09" (r = 0.75") covers a physi- 
cal scale of 0.378-0.763 (3.147-6.354) kpc for these AGNs 
with zw 0.276-3.785. Figures 4(a) and 4(b) show the his- 
tograms of the ratio of counts between these two apertures for 
the AGN hosts for the two ACS bands, respectively. We note 
that these ratios provide a firm upper limit to the AGN con- 
tribution since we did not remove stellar emission within the 
small aperture. For the Veog-band (zgso-band), the mean ratio 
of counts of these AGNs is 0.142 (0.139), and 84.2% (87.2%) 
of them have a ratio less than 0.25, which indicate that the 
optical emission is not dominated by the AGN component. 
Note also that these AGN hosts have count ratios similar to the 
parent population of galaxies. (4) According to Figures 4(c) 
and 4(d), there is no apparent correlation between X-ray lu- 
minosity and the ratio of counts for the AGN hosts, which 
supports the conclusion that the optical emission is not dom- 
inated by the AGN component. Generally speaking, galaxy 
surface brightness dims with increasing redshift and the con- 
trast ratio between a nuclear point source and its host galaxy 
increases toward higher redshifts. This could potentially im- 
pact our conclusion (for high-redshift AGNs) that the optical 
emission is not dominated by the AGN component. However, 
as shown in Figs. 4(c, d), this does not appear to be an issue 
with AGNs in Sample A because the ratio of counts does not 
tend to be larger for sources with larger redshifts. 

Figure 5 shows X-ray luminosity as a function of MIPS 
ALA,24Aim for die 157 AGNs in Sample A that have MIPS 
24 yum detections. There is no apparent correlation between 
Lx and XL\24iim for these AGNs, which is consistent with the 
X-ray and 24 /im emission not being produced by the same 
process, suggesting that, typically, these AGNs do not sig- 
nificantly affect the observed 24 /im fluxes that are used to 
compute the SFRs. 

4.6.4. Sample Overview 

Table 5 presents relevant data and Figure 6 shows his- 
tograms of several physical properties for the sources in 
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Fig. 2. — Plots of Zphot quality checks for Sample A, with AGNs shown as filled circles/solid histogram and non-AGN galaxies shown as g 
histogram, (a) Zp^ot vs. Zspcc. (b) Normalized histograms of Az/(l+Zspec)- The histograms largely overlap each other because AGNs and non-AGN 
photometric redshifts of similar quality, (c) Az/(1 +Zspec) as a function of i?-band magnitude hir. 
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Fig. 3. — (a) 0.5-8 keV luminosity as a function of redshift for Sample A (see § 4.6). Filled circles indicate AGNs and small grey dots indicate non-AGN 
galaxies with their upper Hmits derived from corresponding sensitivity maps. The spread of a factor of Ri 3 in the upper limits at a given redshift reflects the 
sensitivity variation with off-axis angle, i.e., a smaller off-axis angle corresponds to a higher sensitivity. The two dashed lines show the X-ray luminosity cut 
of 41.9 < log(Lx/(ergs"')) < 43.7. For comparison, AGNs that lie outside the luminosity range of 41.9 < log(Lx/(ergs"')) < 43.7 are shown as grey circles, 
(b) 0.5-8 keV luminosity as a function of M\j for the AGNs in Sample A. (c) 0.5-8 keV luminosity as a function of rest-frame U — V color for the AGNs in 
Sample A. AGNs in (b) and (c) are color-coded such that the blue (green, orange, red) color represents 0<z<l(l<z<2, 2<z<3, 3<z<4). No apparent 
correlation is seen in (b) or (c), demonstrating that the optical emission of these AGN hosts should not be dominated by the active nuclei. (A color version of this 
figure is available in the online journal.) 



Sample A. In Sample A, sa 95% (« 98%) of the sources have 
rest-frame /T-band (/-band) coverage or beyond, which en- 
sures the reliability of our stellar-mass estimates (see § 3.2). 
As seen in Fig. 6(c), it is apparent that AGN hosts are gener- 
ally more massive than non-AGN galaxies. 

In Sample A, w 23% (sa 16%) of the galaxies (AGNs) have 
no MIPS 24 fim detection and were given upper limits on 
SFR. At < z < 1 none of the sources without MIPS 24 ^m 
detections has an upper limit on SFR greater than 10 M© yr"', 
and at 1 < z < 2 (2 < z < 3) only two of the galaxies without 
MIPS 24 /im detections have upper limits on SFR greater than 
20 (30) Mq yr"' ; these indicate that Sample A is roughly com- 
plete above SFR« 10 (20, 30) M© yr"' at0<z<l(l<z<2, 
2<z<3). 

Figure 3(a) shows 0.5-8 keV luminosities of AGNs in 
Sample A as a function of redshift. According to Fig. 3(a), 
AGNs in Sample A are complete up to z ~ 1.5, whereas the 
completeness drops toward higher redshifts; above log(Lx) ~ 
43.0, AGNs in Sample A are roughly complete out to z ~ 3. 
We note that the incompleteness of X-ray AGNs at high red- 
shifts does not affect AGN fraction calculations using the 



method detailed in § 5.2.1. 

5. RESULTS 

We chose to use the rest-frame U-V color (i.e., My-My) 
throughout this section because this color straddles the 

o 

4000 A break of the continuum and better tracks the age and 
metallicity variations of the stellar populations in galaxies 
than the U-B (i.e., Mu-Mg) color (e.g.. Bell et al. 2004). We 
have confirmed that the same basic results can be obtained by 
using the rest-frame U-B color 

5.1. Color Bimodality 

In this sub-section, we examine whether there is color bi- 
modality (i.e., separation of the red sequence and the blue 
cloud) for non-AGN galaxies and AGN hosts in Sample A. 

5.1.1. Color Bimodality: Non-AGN Galaxies 

Figure 7 shows histograms of rest-frame U-V color for 
non-AGN galaxies in Sample A, which spans a broad range 
of about -1.0 <U-V < 2.0. The histograms were made with 
a set of redshift bins of Az = 1 .0 over z = 0^. The choice of 
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Table 5 
Relevant Data for the Sources in Sample A 



No. RAj2ooo.o DECj2ooo.o 
(1) (2) (3) 



(4) 



niR 


Mu 


Mb 


My 


log(MB_/v/R/i/j/H/K.* 7^0 ) 


(5) 


(6) 


(7) 


(8) 


(9) 



Fmass SFR(M0 yr-i) XID log[Lx/(ergs-i)] 
(10) (11) (12) (13) 



1 189.37823 62.18428 

2 189.37852 62.18563 

3 189.38080 62.18413 

4 189.26192 62.18623 

5 189.26080 62.18484 



2.486 25.33 -21.73 -21.68 

0.681 22.13 -20.30 -20.09 

1.169 24.49 -20.38 -20.39 

+0.559 22.20 -19.77 -20.01 

+1.014 23.72 -20.50 -19.94 



-22.01 (10.12,10.13,10.11,10.08,9.95,9.89,9.91) 7 

-20.53 (9.74,9.75,9.76,9.80,9.71,9.69,9.63) 7 

-20.80 (9.78, 9.79, 9.77, 9.74, 9.63, 9.59, 9.61) 7 

-20.80 (10.49, 10.50, 10.53, 10.56, 10.46, 10.43, 10.39) 7 

-20.42 (9.76, 9.76, 9.79, 9.85, 9.77, 9.75, 9.70) 7 



8.7 





-42.85 


-1.8 





-41.46 


6.0 





-42.02 


6.0 





-41.03 


3.5 





-41.67 



This table is available in its entirety in a machine-readable form in the online journal. A portion is shown here for guidance regarding its form and content. The 
full table contains 4941 entries. Columns: (1) Source sequence number. (2, 3) J2000 right ascension and declination of the optical counterpart (in degrees), 
respectively. (4) Spectroscopic redshift (with a plus sign) or photometric redshift (without a plus sign). (5) Apparent S-band magnitude (AB mags). (6, 
7, 8) Rest-frame absolute Johnson U-, B-, and V-band magnitude (Vega mags), respectively. (9) Stellar-mass estimate derived at the B-, V-, R-, /-, J-, H-, 
and X'-band (irregardless of whether there is corresponding rest-frame coverage of a source), respectively. For the 4941 sources in Sample A, the median of 
(Mb;v/r/i/]/h.*)/Mk.* is ~ (1.46, 1.49, 1.48, 1.51, 1.19, 1.08), respectively. (10) Flag of reddest rest-frame coverage of a source. This flag shows which 
stellar-mass estimate is adopted as the final estimate for a source (see § 3.2). Fmass =(1, 2, 3, 4, 5, 6, 7) means that Mb/v/r/i/j/h/k.* is adopted for a source, 
respectively. (11) Star-formation rate (positive value) or upper limit on star-formation rate (negative value). (12) X-ray ID from Alexander et al. (2003) 
for the CDF-N AGNs (XID> 0) and Luo et al. (2008) for the CDF-S AGNs (XID< 0), respectively (XID=0 for non-AGN galaxies). Four CDF-N AGNs 
have XID> 503 and are from the CDF-N supplementary optically bright Chandra source catalog; their coiresponding XIDs in this supplementary catalog are 
XID-503. Two CDF-S AGNs have 462 <XID< 548 (XID> 548) and are from the supplementary CDF-S plus E-CDF-S Chandra source catalog (the CDF-S 
supplementary optically bright Chandra source catalog); their corresponding XIDs in this supplementary catalog are XID-462 (XID-548). (13) Intrinsic 
0.5-8 keV X-ray luminosity (positive value; for AGNs) or upper limit on intrinsic 0.5-8 keV X-ray luminosity (negative value; for non-AGN galaxies). 
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Ratio cj V band counts Patio ot z bond counts 

Fig. 4. — AGN contribution to the total (galaxy+AGN) optical emission in 
the Vgoe-band and zgsg-band (see § 4.5) for Sample A. The ratio of counts is 
computed by measuring optical counts in circular apertures of two different 
sizes (r = 3 and 25 pixels, i.e., r = 0.09" and 0.75". respectively) located at 
the centroid of the optical emission. These ratios should provide a firm up- 
per limit to the AGN contribution since we did not remove stellar emission 
within the small aperture, (a, b) Histograms of the ratio of counts of the AGN 
hosts (solid histograms) and non-AGN galaxies (dotted histograms), where 
the non-AGN galaxy numbers have been scaled to match those of the AGN 
hosts. The majority of the AGN hosts have > 75% of their total optical emis- 
sion (in the Vjoe and zgso bands) from outside of r = 0.09" and have count 
ratios similar to the parent population of galaxies, (c, d) X-ray luminosity 
vs. ratio of counts for the AGN hosts in the Vjoe-band and zgso-band. There 
is no clear correlation between X-ray luminosity and the ratio of counts for 
the AGN hosts, which supports the conclusion that the optical emission is 
not dominated by the AGN component. We note that there are a few AGN 
hosts that have a ratio of counts greater than 0.5, which should not be a prob- 
lem because their optical images appear normal and the ratio of counts only 
represents an upper limit of AGN contribution to the total optical emission. 
AGNs in (c) and (d) are color-coded in the same way as those in Figs. 3(b, c). 
(A color version of this figure is available in the online journal.) 



Az =1.0 facilitates examination of possible color evolution 
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Fig. 5. — X-ray luminosity as a function of MIPS AL;^ 24fjm for 'hs 157 
AGNs in Sample A that have MIPS 24 /tm detections. There is no apparent 
correlation between Lx and AL;^ 24fmi. which suggests that these AGNs do 
not significantly affect the observed 24/im fluxes that are used to compute 
the SFRs. 



over cosmic time, avoids likely dilution of potential bimodal 
color behaviors when a wide range of redshift is considered 
(e.g., A2 = 2.0), and avoids poor statistics due to low source 
counts in a narrow redshift bin (e.g., Az = 0.2). 

We find that, in all four redshift ranges, the histogram of 
U -V color for non-AGN galaxies can be acceptably fitted 
(using x^ fitting) by two Gaussian components each in the 
form of A^ = Gn exp(-[(X - Gc)lG^Y/2), where X = U-V, 
Gn is a normalization factor, and Gc (Gw) is the centroid 
(width) of the Gaussian component. ^^ The best-fit Gaussian 
parameters can be found in Table 6. We note that a bimodal 
signature of the galaxy color distribution seems to be present 
in all four redshift ranges, although only a weak redward tail 

^' Based on previous works that established galaxy color bimodality, we 
applied an additional constraint of 0.1 < Gw < 0.5 to the Gaussian fits, i.e., 
we required that the width of the Gaussian component not be too small or 
too large. We did not apply this additional constraint when examining pos- 
sible color bimodality for AGN hosts because there is no established color 
bimodality for AGN hosts. 
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Fig. 6. — Histograms of physical properties for Sample A, with AGNs 
shown as solid lines and non-AGN galaxies as dotted lines: (a) Redshift; (b) 
S-band magnitude; (c) Stellar mass; (d) 0.5-8 keV X-ray luminosity; and (e) 
Effective photon index. The number of non-AGN galaxies has been rescaled 
to match that of AGNs in each plot. The histogram of SFR is not shown, due 
to the existence of upper limits on SFR (see § 3.3). 



can be seen in the highest range (see Fig. 7). We then tested 
for bimodality of the galaxy color distribution using the KMM 
mixture-modeling algorithm (e.g., Ashman, Bird, & Zepf 
1994). This algorithm indicates bimodal color distributions 
with high statistical significance in all four redshift ranges 
shown in Fig. 7, with the derived KMM P-values all being 
equal to zero (see Table 6 for the P-values and corresponding 
mean values of two potential sub-populations). However, we 
caution that the above seemingly established color bimodal- 
ity is based on an assumption of Gaussian-distributed popula- 
tions, which is unlikely to be the true distribution. 

To assess further the validity of the claim of color bimodal- 
ity, we imposed three additional necessary conditions (i.e., 
these three conditions must be satisfied in o rder for a bi modal- 

ity claim): (1) The quantity |Gci-Gc2|/a/ Gwi + ^W2 should 

be equal or greater than 0.707 (i.e., -\/2/2), where the value 
of 0.707 is obtained when |Gci -Gcil = Gw\ = Gwi- This 
quantity is a measure of separation between the centroids of 
two components. (2) The quantity Gni/Gn2 should be in a 
range roughly from 0.1 to 10, i.e., one component should not 
be negligible compared to the other. (3) The requirement of 
a two-Gaussian-component fit (rather than a single-Gaussian- 
component fit) is statistically significant. To address Condi- 
tion (3), we calculated the F-statistic and its associated proba- 
bility given "new" and "old" values of x^ and degrees of free- 
dom (i.e., values derived from two-Gaussian-component fits 
vs. those derived from single-Gaussian-component fits). A 
low probability from the F-test indicates that two-Gaussian- 
component fits are better than single-Gaussian-component fits 
at a high confidence level. As shown in Table 6, the cases of 



< z < 1, 1 < z <2, and 2 < z < 3 that have apparent bi- 
modal signatures (see Fig. 7) satisfy the above three condi- 
tions, while the case of 3 < z < 4 that has a less prominent 
bimodal signature (also see Fig. 7) fails the third condition 
(i.e., the requirement of a two-Gaussian-component fit is not 
statistically significant). 

As a final check, we calculated the separation lines be- 
tween the red sequence and the blue cloud for different red- 
shift ranges using the following equation and plotted them as 
vertical solid lines in Fig. 7: 



(f/ - V)re,st = -0.3 lz-0.08Mv- 0.51. 



(4) 



This equation was derived by Bell et al. (2004) who stud- 
ied the color distribution of «25,000 R <24 galaxies with 
0.2 < z < 1.1. These four separation lines shown in Fig. 7 
(median values of z and My used for each range of redshift) 
seem to describe our observational data reasonably well, tak- 
ing into account a typical color scatter of < 0.2 mag for the 
red sequence color-magnitude relation (see § 4 in Bell et al. 
2004). The general agreement between the separation lines 
defined by Eqn. (4) and our data suggests: (1) the term of 
-O.OSMv in Eqn. (4) may be valid up to at least z « 3 since 
the range of My of the galaxies in Sample A is similar to that 
of the Bell et al. galaxies; (2) the term of -0.31z may also be 
valid up to at least z ~ 3, which indicates that the definitions 
of the red sequence and the blue cloud may evolve with red- 
shift, i.e., galaxies may tend to be bluer as they are younger 
(i.e., at higher redshifts; also see § 5.3). We note that a care- 
ful extension of Eqn. (4) out to z ~ 4 would be useful, but is 
beyond the scope of this work. 

Based on the above analyses, we conclude that there is ev- 
idence that galaxy color bimodality holds from the local uni- 
verse up to z ~ 3, with likely hints of color bimodality up to 
z « 4. We note that the relative fraction of galaxies in the 
red and blue populations evolves strongly, i.e., the fraction of 
red-sequence galaxies generally decreases as the redshift in- 
creases (we caution that this trend is affected quantitatively 
by the incompleteness of red galaxies in Sample A at high 
redshifts); and that the colors of both the blue and red com- 
ponents become redder from the early universe to the present 
time (see Fig. 7 and the values of Gni/Gn2, Gqi, and Gc2 in 
Table 6). Our result confirms many previous works on galaxy 
color bimodality (e.g., Strateva et al. 2001, Hogg et al. 2002, 
and Baldry et al. 2004 for the local universe; Bell et al. 2004, 
Weiner et al. 2005, Willmer et al. 2006, Cirasuolo et al. 2007, 
Franzetti et al. 2007, and Taylor et al. 2009a for up to z ~ 1-2; 
Kriek et al. 2008 for the detection of a red sequence of mas- 
sive field galaxies at z ~ 2.3; Brammer et al. 2009 for up to 
z ~ 2.5; Giallongo et al. 2005 for up to z ~ 2.5-3) and extends 
them to higher redshifts. However, we note that some authors 
found no color bimodality for galaxies at high redshifts (e.g., 
Cirasuolo et al. 2007 for z > 1.5; Labbe et al. 2007 for no 
well-defined red sequence at z ~ 3). Therefore, the situation 
is controversial at z > 1.5, unlike the well-established galaxy 
color-bimodality below z ~ 1.5. 

5.1.2. Color Bimodality: AGN Hosts 

We performed the same procedure as for the case of non- 
AGN galaxies to test for color bimodality of AGN hosts in 
Sample A. The results are shown in Fig. 8 and Table 6. We 
note that only two redshift ranges were considered (0.0 < z < 
1.0 and 1.0 < z < 2.0); the redshift ranges of 2.0 < z < 3.0 
and 3.0 < z < 4.0 were not considered due to an insufficient 
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-1.0 0.0 

U—V (rest — frame) 

Fig. 7. — Histograms of rest-frame U — V color for non-AGN galaxies in Sample A (shown as black solid lines) in different redshift intervals, with a bin size 
of A(t/-y) = 0.1. The peak value of each histogram (Wpeak) has been rescaled to unity and the total number (Mot) of non-AGN galaxies in each subsample is 
shown. Thick blue dotted lines show the two best-fit Gaussian components, with their sum shown as the red solid line. Best-fit Gaussian parameters and relevant 
statistical properties are shown in Table 6. A set of separation lines between the red sequence (to the right) and the blue cloud (to the left) is shown as vertical 
solid hues that are derived using Eqn. (4) (see § 5.1.1). As a comparison, a set of lines representing the color-magnitude relation of red-sequence galaxies is 
shown as vertical dashed lines, which are roughly 0.25 mag redward of the separation lines (Bell et al. 2004). (A color version of this figure is available in the 
online journal.) 



Table 6 
Best-Fit Gaussian Parameters and Relevant Statistical Properties for Figs. 7, 8, & 9 









Gaussian 1 


Gaussian! KMM-test |Cci-Gc2|/ /^-test 


Figure 


z 


Mot Mpeak 


(Gci, Cwi, C?Ni) 


(Gc2, Gw2, Gn2) (/'-value, mean 1, mean 2) (G^j +G|,2)°'^ Gni/Gn2 P 


(1) 


(2) 


(3) (4) 


(5) 


(6) (7) (8) (9) (10) 



Fig. 7(a) 0.0<z< 1.0 1626 

Fig. 7(b) 1.0<z<2.0 1570 

Fig. 7(c) 2.0<z<3.0 1161 

Fig. 7(d) 3.0<z<4.0 584 



230 (0.490,0.213,0.426) (1.088,0.209,0.716) 

195 (0.242,0.231,0.753) (0.955,0.248,0.439) 

200 (0.052,0.251,0.762) (0.824,0.196,0.094) 

133 (-0.156,0.178,0.667) (0.465,0.250,0.079) 



(0.000,0.478,1.112) 2.01 0.60 0.0381 

(0.000,0.272,0.981) 2.10 1.72 0.0029 

(0.000,0.045,0.799) 2.42 8.14 0.0256 

(0.000,-0.189,0.600) 2.02 8.45 0.2407 



Fig. 8(a) 0.0<z< 1.0 83 17 
Fig. 8(b) 1.0<z<2.0 85 11 



(1.014, 0.481, 0.318) (1.091, 0.056, 0.899) 
(0.537, 0.554, 0.259) (0.982, 0.205, 0.681) 



(0.003,0.489,1.081) 0.16 0.35 0.0513 

(0.002,0.130,0.940) 0.75 0.38 0.2866 



Fig. 9(a) 0.0<z< 1.0 1626 

Fig. 9(b) 1.0<z<2.0 1570 

Fig. 9(c) 2.0<z<3.0 1161 

Fig. 9(d) 3.0<z<4.0 584 

Fig. 9(A) 0.0<z< 1.0 83 

Fig. 9(B) 1.0<z<2.0 85 



197 (0.088,0.167,0.561) (0.759,0.273,0.826) 

179 (0.066, 0.204, 1.001) (0.727, 0.266, 0.510) 

232 (-0.042, 0.225, 0.733) (0.648, 0.214, 0.096) 

125 (-0.244, 0.222, 0.695) (0.466, 0.250, 0.039) 

14 (0.069, 0.250, 0. 152) (0.790, 0.237, 0.852) 

16 (0.450,0.500,0.360) (0.751,0.052,0.700) 



(0.000,0.153,0.804) 2.09 0.68 0.0000 

(0.000,0.099,0.778) 1.98 1.96 0.0000 

(0.000,-0.042,0.699) 2.22 7.61 0.0061 

(0.000,-0.252,0.600) 2.13 17.85 0.3220 

(0.002,0.056,0.789) 2.09 0.178 0.0385 

(0.003,0.470,0.553) 0.60 0.51 0.0117 



Columns: (1) Subplot considered. (2) Redshift range considered. (3) Total source number in each subplot. (4) Peak value of the histogram in each subplot. (5, 6) 
Gaussian component (centroid, width, normalization) that corresponds to the blue and red population, respectively. (7) Result of the KMM mixture-modeling 
algorithm that tests for bimodality. A P-value of < 0.05 is required for a bimodafity claim. Mean 1 and mean 2 are the mean values of two potential 
sub-populations. (8) Quantity that measures the separation between the centroids of two Gaussian components. A value of < 0.707 is required for a bimodality 
claim. (9) Quantity that measures the relative strength between two Gaussian components. A value ranging from 0.1 to 10 is required for a bimodality claim. 
(10) Probability from the F-statistic that measures whether the requirement of a two-Gaussian-component fit rather than a single-Gaussian-component fit is 
statistically significant. A value of < 0.05 is required for a bimodality claim. 



number of sources for meaningful analysis. For 0.0 < z < 1.0 
the condition of |Gci 



-Gc2|/ 



G^[ -I- G\i2 ^ 0.707 is not sat- 
isfied (a value of 0.16 was obtained) and the requirement of 
a two-Gaussian-component fit rather than a single-Gaussian- 
component fit is not statistically significant (p = 0.0513 in the 
F-test). For 1 .0 < z < 2.0, the requirement of a two-Gaussian- 
component fit is not statistically significant {p = 0.2866 in the 
F-test). Therefore, our result reveals no apparent color bi- 
modality for AGN hosts, which is in agreement with previ- 
ous studies (e.g., Bohm & Wisotzki 2007 for a sample with 
a mean redshift of (z) « 0.6; Silverman et al. 2008b for 
0.4 < z < 1.1) and extends them to z ~ 1-2. 

According to Fig. 8, it also seems clear that AGN hosts are 
generally redder than non-AGN galaxies up to z ~ 3^. This 
result confirms previous works (e.g., Barger et al. 2003; Nan- 



dra et al. 2007; Brusa et al. 2009; Treister et al. 2009; but see 
§ 5.3 and 5.4). 

5.1.3. Color Bimodality: Effects of Dust Reddening 

Recent studies have investigated the effects of dust red- 
dening when analyzing galaxy colors (e.g., Cowie & Barger 
2008; Brammer et al. 2009) and find that many of the green- 
valley sources in the CMD are dust-reddened blue-cloud 
sources, rather than sources undergoing a blue-to-red transi- 
tion. Specifically, Brammer et al. (2009) studied a /T-band 
selected (K < 22.8) galaxy sample and argued that correcting 
the rest-frame colors for dust reddening allows an improved 
separation between the red sequence and the blue cloud up 
to z ~ 2.5. To examine the effects of dust reddening on 
Sample A, we plot the dust-corrected color (i.e., L'^-V- Auv) 
histograms for both non-AGN galaxies and AGN hosts in 
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-1.0 0.0 

U—V (rest — frame) 

Fig. 8. — Same as Fig. 7, but for AGN hosts (shown as black sohd hnes) in Sample A. For comparison, histograms for non-AGN galaxies (from Fig. 7) ai'e also 
shown as the dotted lines. Note that Gaussian fits were applied only to Panels (a, b) which have a sufficient number of AGNs for meaningful analysis. Best-fit 
Gaussian parameters and relevant statistical properties for these two panels are also shown in Table 6. AGN hosts are generally redder than non-AGN galaxies. 
(A color version of this figure is available in the online journal.) 




0.0 1.0 -1.0 0.0 1.0 

t/— K— Ayy (rest — frame) 

Fig. 9. — Same as Figs. 7 and 8, but correcting for dust extinction for non-AGN galaxies (Panels (a - d)) and AGN hosts (Panels (A - D)) in Sample A, 
respectively. Best-fit Gaussian parameters and relevant statistical properties ai'e also shown in Table 6. (A color version of this figure is available in the online 
journal.) 



Sample A in Fig. 9. For this assessment, we followed Eqn. (2) 
of Brammer et al. (2009), Auv = 0.47Av, to correct for dust 
reddening. We obtained V-band extinction information (i.e.. 
Ay) from the template SED fitting (see § 3.1). 

According to Fig. 9 and Table 6, it seems, after correct- 
ing for dust extinction, that: (1) For non-AGN galaxies at 
0<z<l,l<z<2, and 2 < z < 3, there is apparent improve- 
ment of the separation between the red sequence and the blue 
cloud, in the sense that the profiles of the red and blue compo- 
nents appear smoother and clearer and the p values from the 
F-test become significantly smaller; but for non-AGN galax- 
ies at 3 < z < 4, there is still no apparent color-bimodality 
(both the value of Gni/Gn2 and the p value from the F-test 
are too large). (2) For AGN hosts at < z < 1, there appears 
to be some improvement of the separation between the red 



sequence and the blue cloud. Following the procedure de- 
scribed in § 5.1.1, we find likely evidence of color bimodality 
for AGN hosts at < z < 1 (we caution that the relatively 
weak blue component could simply be a blueward skewed 
tail of the red component). For AGN hosts at 1 < z < 2, we 
still find no evidence of color bimodality (the condition of 

|Gci -Gci]/ \/G^i + G^2 — 0.707 is not satisfied; a value of 
0.60 was obtained). (3) Generally, AGN hosts still appear red- 
der than non-AGN galaxies at all redshift ranges (comparing 
Figs. 9(a - d) with Figs. 9(A - D), respectively). 

We find that there are systematic blueward shifts in the col- 
ors due to de-reddening by comparing Figs. 7 and 8 with 
Fig. 9, which suggests that dust extinction could be responsi- 
ble for the red colors of dusty star-forming galaxies including 
many AGN hosts (also see, e.g., Brusa et al. 2009). How- 
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ever, based on the above analysis, we conclude that dust ex- 
tinction should not affect our results on color-bimodality in 
a significant way. We have also verified that the results ob- 
tained hereafter will be unchanged qualitatively after using 
the above dust de-reddening technique. Furthermore, most 
previous CMD-related works are based on dust-uncorrected 
colors, e.g., the separation line between the red sequence and 
the blue cloud (i.e., Eqn. (4)). Therefore, we adopted dust- 
uncorrected colors hereafter for the purpose of straightfor- 
ward comparison with previous works. We note that a more 
careful assessment of the effects of dust reddening is desirable 
in the future, but is beyond the scope of this work. 

5.2. Color-Magnitude Relations 

In this section, we first investigate color-magnitude rela- 
tions of AGN hosts and non-AGN galaxies; we then explore 
the location of luminous submillimeter galaxies in the CMD 
and color-mass diagram. 

5.2. 1 . CMDs for AGN Hosts and Non-AGN Galaxies 

The CMDs for the AGNs and their parent galaxies in Sam- 
ples B, C, and D are shown in Figs. 10(a - i). In each CMD, 
we divided sources (both non-AGN galaxies and AGN hosts) 
into three groups (the red sequence, the green valley, and the 
blue cloud) using Eqn. (4) on a source-by-source basis. If a 
source satisfies (f/-y),.est+0.31z+0.08Mv+0.51 > 0.05, it is 
in the red sequence; if a source satisfies -0.05 < (f/- V)rest + 
0.3 lz+0.08Mv + 0.51 < 0.05, it is in the green valley; and if 
a source satisfies (f/ - y)rest + 0.3 lz+0.08Mv + 0.51 < -0.05, 
it is in the blue cloud. Figures 10(A - I) show corresponding 
plots of AGN fraction as a function of rest-frame U-V color 
for Samples B, C, and D. For a given X-ray luminosity range, 
we computed the AGN fraction (/) and its associated error 
(ct/) according to the following two equations that were used 
by Silverman et al. (2008): 



^ 1 



-}" 



E 






(5) 



(6) 



where / and a/ are a sum over the full sample of AGNs (A^) 
with A^gai.i denoting the number of galaxies capable of host- 
ing the /th detectable AGN with X-ray luminosity L^- This 
method not only takes into account the spatially varying sen- 
sitivity limits of the 2 Ms CDF observations using the sen- 
sitivity maps derived in § 4 of Luo et al. (2008), but also 
treats the incompleteness of X-ray AGNs at high redshifts [see 
Fig. 3(a)] effectively so that it does not affect AGN fractions. 
We have the following observations up to z w 2-3 of the 
material in Fig. 10. (1) The majority of AGNs reside in 
massive galaxies (comparing AGN counts in Figs. 10(a - c) 
with those in Figs. 10(d - f); also see Fig. 1 and Fig. 6(c); 
see § 5.3 for further relevant details). (2) For Sample C 
(i.e., M* > 1O'"^M0), AGN hosts do not appear more lu- 
minous than non-AGN galaxies; both AGN hosts and non- 
AGN galaxies are luminous and share similar My distribu- 
tions. This is in contrast to the result from either Sample B 
(i.e., M* > IO'^-^Mq) or Sample D (i.e.. My < -19) that AGN 
hosts are generally more luminous than non-AGN galaxies 
(i.e., having more negative values of My); the latter result was 
also found by past CMD works (e.g., Nandra et al. 2007; 



Silverman et al. 2008b). (3) For Sample C, AGNs do not 
seem to reside predominantly in the red sequence, the top of 
the blue cloud, or the green valley in between; in fact, AGN 
hosts seem to spread throughout the CMDs (see, particularly. 
Figs. 10(d, e)). This is different from the result obtained with 
either Sample B or Sample D that there is apparent AGN clus- 
tering in the CMDs; the latter result was also obtained by 
previous CMD studies (e.g., Nandra et al. 2007; Rovilos & 
Georgantopoulos 2007; Westoby et al. 2007; Silverman et al. 
2008b; Hickox et al. 2009; Schawinski et al. 2010). (4) For 
Sample C, the AGN fraction seems to remain broadly con- 
stant (at « 10%) regardless of color [see Figs. 10(D - F)]. 
This is in contrary to the case of either Sample B or Sample D 
where the AGN fraction generally increases as colors become 
redder, with likely signs of leveling off at red colors [see 
Figs. 10(A - C, G - I)]. Whether a large population of low- 
mass blue galaxies is included in the analysis holds the key 
to explaining the above different trends of AGN fraction seen 
between Sample C and Sample B/D, i.e., this population of 
galaxies is included in Sample B/D, which is responsible for 
the drop in AGN fractions at blue colors (also see Silverman 
et al. 2009). Results on the fraction of X-ray luminous AGNs 
(i.e., Lx > lO'*-' erg s"') as a function of color (for Samples B, 
C, and D) are in general agreement with the above results on 
the AGN fraction. We note that the results on AGN fraction 
for Sample C likely indicate the duty cycle of SMBH growth 
in typical massive galaxies. (5) For Samples B, C, and D, it 
seems that the majority of X-ray hard AGNs (i.e., Feff < 0.5) 
reside in the red sequence, which we suspect is due to the fact 
that the majority of AGNs overall reside in the red sequence. 
Indeed, Fisher's exact probability test shows no strong evi- 
dence [values of Ppisher > 0. 1 are obtained except for one case 
(Fig. 10(a)) where Ppisher = 0.034] that X-ray hard AGNs pref- 
erentially reside in the red sequence compared to AGNs over- 
all. We note that the incidence of X-ray hard AGNs in the red 
sequence was reported by previous works and has been used 
to argue for BH accretion persisting after the termination of 
star formation (e.g., Nandra et al. 2007; Georgakakis et al. 
2008). 

5.2.2. CMDs for Luminous Submillimeter Galaxies 

Luminous submillimeter galaxies (SMGs) are a class of 
strongly star-forming, dust-obscured galaxies at z « 1-3 that 
have been studied extensively (e.g., Alexander et al. 2005a, 
2005b; Chapman et al. 2005; Pope et al. 2006, 2008; Biggs 
et al. 2010; Laird et al. 2010; Swinbank et al. 2010; Wardlow 
et al. 2010) but have yet to be put into the context of CMDs. 
Since we can now extend CMD studies of AGN hosts and 
non-AGN galaxies into this redshift regime, it is of interest to 
see where SMGs lie in CMDs, which may provide one useful 
guide to interpreting CMD results. 

The accuracy of submillimeter positions (for SMGs) is rela- 
tively poor (typically accurate to « 1-7"; see, e.g.. Biggs etal. 
2010), compared to those of optical, infrared, or X-ray posi- 
tions. Thus it is challenging to identify secure counterparts for 
SMGs. Chapman et al. (2005) and Biggs et al. (2010) have 
identified reliable radio and/or MIPS 24 /im counterparts for 
samples of SMGs in the CDF-N and E-CDF-S (i.e., Extended- 
Chandra Deep Field-South; Lehmer et al. 2005), respectively, 
through careful analyses. We therefore cross-matched the 
optical positions of the sources in Sample A with the radio 
and/or MIPS 24 ^m positions of the CDF-N and E-CDF-S 
SMGs provided by Chapman et al. (2005) and Biggs et al. 
(2010), respectively (a maximum matching radius of 0.75" is 
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Fig. 10. — (a — i) Color-magnitude diagrams for Samples B. C, and D: rest-frame U -V colors are shown against V-band absolute magnitudes, with AGN 
hosts shown as filled circles and non-AGN galaxies as small dots. Large filled circles indicate X-ray luminous AGNs with Lx > 10*' erg s"' . Diamonds indicate 
X-ray hai'd AGNs with F^ff < 0.5. Red/green/blue symbols represent sources in the regions of red sequence/green valley/blue cloud (see § 5.2). Note that the 
overlaps between the regions of red sequence/green valley/blue cloud are caused by our source-by- source classification scheme (see § 5.2). The solid lines show 
the color limits that separate the red sequence and the blue cloud, derived using Eqn. (4) and the median redshift of galaxies for each subsample. The dashed 
lines correspond roughly to the limit imposed by the stellar-mass cut of either M* > IO'^Mq or M* > 10'" -'Mq. The dotted lines roughly show the flux limit 
of mR < 26. Note that the apparent horizontal stripes in the plots are due to the fact that the galaxies on the same stripe have the same (or very similar) best-fit 
template but different redshifts, which leads to the same (or very similar) U -V colors but different values of M\ . (A - 1) AGN fraction as a function of rest-frame 
U — V color (in bins of A(J/ — V) = 0.4) corresponding to each of the above cases of (a - i), respectively. Black solid (dashed-dot) curves show the fractions of 
AGNs (X-ray luminous AGNs; i.e., 43.0 < log(Lo5-8kcV /(£■■§ s"')) < 43.7); orange (purple) filled patterns show the corresponding \cr confidence ranges for 
AGNs (X-ray luminous AGNs); and orange (purple) numbers show the ranges of X-ray luminosity for AGNs (X-ray luminous AGNs) considered. The vertical 
lines roughly show the separation lines derived using Eqn. (4) and the median redshift and M\ of galaxies for each case. (A color version of this figure is available 
in the online journal.) 
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adopted; the associated false-match probability is « 3.9%). 
We identified 11 SMGs with 1.0 < z < 3.0 (with a median 
redshift of « 2) in Sample A; 5 of these SMGs host AGNs.^*" 
Figure 1 1 shows the location of these 1 1 SMGs in the CMD 
and color-mass diagram. These SMGs are spread throughout 
the CMD in the sense that they are found in the red sequence, 
the blue cloud, and the green valley in between, although on 
average they appear more luminous than the general galaxy 
population (SMGs with SFRs of « 500 Mq yr"' would have 
Luv ~ 5 X IO'^Lq if they were not obscured). This is likely 
due to extinction effects in these extreme systems: the red 
colors might well indicate dust extinction rather than old stel- 
lar populations; and the blue colors might reflect that a small 
fraction of the star formation is unobscured and thus domi- 
nates over the obscured star-formation component in the ob- 
served optical-near-lR band. Therefore, the observed colors 
of these SMGs do not represent their intrinsic colors. In the 
color-mass diagram, these SMGs, as expected, lie toward the 
high stellar-mass end (i.e., they are much more massive than 
the general galaxy population); the median stellar mass of 
these SMGs is « 1.6 x IO'^Mq, which is toward the low end 
of the estimates from Borys et al. (2005). The location of 
SMGs in the CMD and color-mass diagram provides a pow- 
erful and specific example of how CMD and color-mass di- 
agram results can be subject to reddening effects, which un- 
derscores the need for caution when interpreting the color- 
magnitude properties of small samples of extreme systems. 

5.3. Color Dependence on Physical Properties 

Measurements of the dependences of galaxy colors on 
physical properties such as redshift, SFR, and stellar mass 
can provide constraints on models of galaxy formation and 
evolution. In this sub-section, we examine the dependences 
of both AGN host colors and non-AGN galaxy colors on red- 
shift, SFR, and stellar mass, and also examine the likely cor- 
relations between them. The dependences of AGN host col- 
ors and non-AGN galaxy colors on redshift, SFR, and stellar 
mass for Sample B are shown in Fig. 12 and the correspond- 
ing statistical properties are shown in Table 7. Based on the 
information in Fig. 12 and Table 7, we obtain the following 
results: 

(1) AGN hosts are generally redder than non-AGN galaxies 
irrespective of redshift or SFR range considered [see the run- 
ning medians of colors plotted in Fig. 12(a), the color distribu- 
tions in Figs. 12(b - d, f, g), and the Pks values^^ in Table 7]. 
However, we find that this trend seems to be valid only up to 
M* K, 1O'*'^M0 when examining the colors as a function of 
steUar mass; above IO^'-'Mq, both AGN hosts and non-AGN 
galaxies have similar color distributions [see the running me- 
dians of colors plotted in Fig. 12(h), the color distributions 
in Figs. 12(i - k), and the Pks values in Table 7]. This lat- 
ter observation appears to be inconsistent with the result from 

^* For the CDF-N, Chapman et al. (2005) presented a spectroscopically 
identified sample of 22 SMGs, among which 19 are located within the central 
''encircled =8' radius area; of these 1 9 SMGs, 1 5 have counterparts in our North 
base catalog and only 11 are in Sample A. For the E-CDF-S, Biggs et al. 
(2010) presented a submillimeter-flux-limited sample of 126 SMGs, among 
which 12 are located within the central rcncircled =8' radius area; of these 
12 SMGs, 5 have counterparts in our South base catalog and only 2 are in 
Sample A. We thus have a total of 13 SMGs in Sample A; 1 1 of them have 
1.0<z<3.0. 

^^ Pks is the probability that the two color distributions considered could 
be drawn from the same parent population, according to the null hypothesis 
from a Kolmogorov-Smimov (K-S) test; a small value of Pks (eg-, < 5.0%) 
indicates that the two color distributions considered are different. 



§ 5.1 that non-AGN galaxies are bimodal in color but AGN 
hosts are not (see Fig. 7 and Fig. 8). We note that this appar- 
ent inconsistency is due to the fact that samples with different 
steflar-mass cuts (i.e., M* > IO'^'^Mq vs. M* > lO'^'^M©) 
were used respectively. We thus made versions of Fig. 7 and 
Fig. 8 for Sample C (i.e., with M* > lO'^^M©) and find that 
the color distributions of non-AGN galaxies and AGN hosts 
are similar and neither shows apparent bimodality. 

(2) There is an apparent color evolution as a function of red- 
shift such that the higher the redshift, the bluer the colors, for 
both AGN hosts and non-AGN galaxies (see Figs. 12(a-d) 
and the values of Spearman's rank correlation coefficient (rs) 
as well as the associated probability (/?s)^** in Table 7). This 
result does not appear to be a selection effect because it is also 
obtained when examining Sample C that is roughly complete 
for both blue and red galaxies. 

(3) From Fig. 12(e), it seems that the AGN fraction rises to- 
ward higher SFRs (also see Rafferty et al. 2010). To examine 
this result further, we plot AGN fraction as a function of SFR 
in different ranges of redshift for Sample B and Sample C in 
Fig. 13. For Sample B, we find that the fraction of AGNs gen- 
erally rises as the SFR increases (note that at the high-SFR 
end this trend is roughly consistent with the AGN fraction 
for z « 2 SMGs obtained by Alexander et al. 2005b). This 
trend persists if we consider the X-ray luminous AGNs [i.e., 
43.0 < log(LA'/(erg s"')) < 43.7]. This result impHes that host 
galaxies of X-ray AGNs generally have higher SFRs than non- 
AGN galaxies from z ~ 0-3 for Sample B. However, when 
examining Sample C, we find that the above trend of increas- 
ing AGN fraction toward higher SFRs becomes less promi- 
nent at 1 < z < 2 and is not detectable at 2 < z < 3, although 
this trend does exist at < z < 1; similar results were ob- 
tained considering the X-ray luminous AGNs. The decrease 
in the strength of the AGN fraction vs. SFR correlation with 
increasing redshift, at least from z ~ 0-2, appears to be a real 
effect not attributable to limited source statistics; this is dis- 
cussed further in § 5.4. 

(4) The presence of about a dozen AGN hosts (with red- 
shifts up to « 1.3 and a median value of Lx ~ 7 x 10^^ erg s~') 
with very low SFRs (i.e., < I.QMq yr"'; see Fig. 12(e)) im- 
plies that AGN activity may persist up to z ~ 1.3 after the 
process of active star formation is largely over, which agrees 
with previous works (e.g., Nandra et al. 2007; Georgakakis et 
al. 2008). 

(5) There is a strong color evolution as a function of stel- 
lar mass such that the more massive the galaxy, the redder 
its color, for both AGN hosts and non-AGN galaxies (see 
Figs. 12(h - k) and the values of rs in Table 7). Furthermore, 
the AGN fraction seems to rise strongly toward higher stel- 
lar masses (see Fig. 12(h)). Figure 14 substantiates these two 
results, showing the rest-frame U -V color and AGN frac- 
tion as a function of stellar mass in different ranges of red- 
shift for Sample B (similar results were obtained consider- 
ing Sample C). First, both AGN hosts and non-AGN galax- 
ies generally become redder as the stellar mass increases ir- 
regardless of redshift range considered. This result agrees 
with previous works at lower redshifts (e.g., Baldry et al. 
2006). The color histograms for both AGN hosts and non- 
AGN galaxies in Sample B and Sample C (shown in the in- 

^** ps is the two-sided significance of deviation from zero and a small ps 
value indicates that the result of rs is obtained at a high significance level 
(i.e., is statistically meaningful). In this paper, we adopt the following cri- 
teria: |rs| ~ 0.0-0.3 indicates no apparent correlations, jrsl ~ 0.3-0.6 indi- 
cates apparent correlations, and | rs | ~ 0.6-1 .0 indicates strong correlations. 
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Fig. 11. — Submillimeter galaxies (shown as crosses) in the color-magnitude diagram [Panel (a)] and color-mass diagram [Panel (b)] for 
1.0 < z < 3.0 (symbols and colors have the same meaning as those in Fig. 10). (A color version of this figure is available in the online journal.) 



10^^ 

Sample A at 



Table 7 
Statistical Properties of Color Distributions for Fig. 



12 



Figure 


AGNs 


{U-V)agn 
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Galaxies 
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K-S (%) 

agn-gal 
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(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 



a) 


185 (100.0%) 


0.89 


-0.33 (0.000) 


4357 (100.0%) 


0.47 


-0.63 (0.000) 


0.0 


0.0<z<3.0 


b) 


83 (44.8%) 
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— 
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and non-AGN galaxies in each subplot, respectively. (4, 7) Spearman's rank correlation coefficient rs and its associated probability ps for AGNs and non-AGN 
galaxies, respectively (see Footnote 28). (8) Probability (Pks) that the color distiibutions of AGNs and non-AGN galaxies could be drawn from the same parent 
population, derived from a Kolmogorov-Smimov (K-S) test (see Footnote 27). (9) Relevant condition appHed in each subplot. (10-12, 13 - 15) The values of 
rsips) derived when combining AGNs or non-AGN galaxies in two subplots, respectively. Note that we did not derive any rs(ps) for SFR-related cases since 
some sources have only upper limits on SFR. 



sets) strengthen a result obtained earlier in this section: for 
Sample B, AGN hosts generally appear redder than non-AGN 
galaxies; whereas for Sample C, AGN hosts and non-AGN 
galaxies have similar color distributions, no matter which 
range of redshift is considered. Second, the fraction of AGNs 
increases strongly as the stellar mass increases irregardless 
of redshift range considered (also true for X-ray luminous 
AGNs), which agrees with the result obtained in § 5.2 that 
most AGNs reside in massive hosts and extends previous find- 
ings to higher redshifts (e.g., Kauffmann et al. 2003; Best et 
al. 2005; Bundy et al. 2008; Alonso-Herrero et al. 2009; 



Brusa et al. 2009; Silverman et al. 2009). We note that for 
a given lower limit on the AGN luminosity, we can identify 
lower Eddington ratio (L/Ledd) AGN activity in higher-mass 
galaxies than in lower-mass galaxies, which introduces an in- 
herent bias toward the above trend of an increasing AGN frac- 
tion with increasing stellar mass. A careful investigation of 
this bias would be useful, but is beyond the scope of this work. 

5.4. Mass-Selection Effects 

In § 5.2, § 5.3, and Fig. 10, we showed that results about 
color-magnitude relations of AGN hosts and non-AGN galax- 



18 



Xue et al. 




-I . ■ I--.--' I 

0.5 
U-V (rest-frame) 



0.0 0.5 1.0 

U-V (rest-frame) 



0.0 0.5 

U-V (rest-frame) 



Fig. 12. — Dependences of color on redshift (left panels), SFR (middle panels), and stellar mass (right panels) for the AGN hosts (filled circles) andnon-AGN 
galaxies (small grey dots) in Sample B. (a) Plot of rest-frame U -V color vs. redshift. (b - d) Histograms of rest-frame U — V color for the AGN hosts (solid 
lines) and non-AGN galaxies (dotted lines) in three ranges of redshift. (e) Plot of rest-frame U -V color vs. SFR. The leftward arrows indicate upper limits on 
SFR. (f, g) Histograms of rest-frame U -V color for the AGN hosts (solid lines) and non-AGN galaxies (dotted lines) in two ranges of SFR. The y-axis in (f) 
has a range of 0-70. (h) Color-mass diagram, (i - k) Histograms of rest-frame U -V color for the AGNs (solid lines) and non-AGN galaxies (dotted lines) in 
three ranges of stellar mass. In panels (a) and (h), the crosses (pluses) show the running medians of U -V colors for AGN hosts (non-AGN galaxies) in bins of 
Az = 0.4 and Alog(Mvt) = 0.5, respectively. The running medians ofU — V colors were not computed in (e) due to the existence of upper limits on SFR. In (a, e, 
h), the vertical dashed lines indicate the parameter values adopted to divide into subsamples. For the histograms, non-AGN galaxy numbers have been rescaled 
to match those of AGN hosts; vertical solid (dotted) lines indicate median AGN host (galaxy) U -V colors of each subsample. (A color version of this figure is 
available in the online journal.) 



ies change when considering Sample B vs. Sample C, which 
have different stellar-mass cuts. For example, the AGN frac- 
tion rises toward redder colors in Sample B, but the AGN frac- 
tion remains nearly constant regardless of color in Sample C. 
We note that over the wide stellar-mass range that makes up 
Sample B {M^, > IO'^'^Mq) there is significant stellar-mass 
dependent variation in the distribution of U -V colors (see 
Fig. 12(h)). Since the AGN fraction increases with stellar 
mass and the U-V color becomes redder with increasing 
stellar mass, we will artificially find that the AGN fraction 
is larger for redder galaxies. Furthermore, the incompleteness 
for red galaxies at high redshifts in Sample B also contributes 
to this bias because the inclusion of a large population of low- 



mass blue galaxies (due to a relatively low stellar-mass cut) 
makes the AGN fraction become smaller toward blue colors. 
In order to mitigate this bias, we have implemented a mass- 
matching technique. We constructed a mass-matched sam- 
ple as follows: for each AGN host in Sample B, ten unique 
galaxies (i.e., no duplicates) with a similar stellar mass (i.e., 
M*.agn/2 < M^,fiAL < 2Mi agn) in Sample B were randomly 
selected (ignoring all information about their My, U-V, and 
SFR values). Note that Sample C happens to be a roughly 
mass-matched sample (see Fig. 12(h)) that ensures a fair com- 
parison between AGN hosts and non-AGN galaxies. 

The top and middle panels of Fig. 15 show CMDs and AGN 
fractions as a function of color for such a mass-matched sam- 
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Fig. 13. — AGN fraction as a function of SFR in bins of Alog(SFR) = 0.5 for Sample B (top panels) and Sample C (bottom panels), respectively (symbols have 
the same meaning as those in Fig. 10). A cut of SFR > 10 (> 20, > 30) Mq yr"' was applied for the redshift bin of 0.0 < z < 1.0 (1.0 < z < 2.0, 2.0 < z < 3.0) 
in order to remove the uncertainties caused by upper Umits on SFR (see § 4.6.4). The filled circles indicate the approximate AGN fraction for z~1 submillimeter 
galaxies (Alexander et al. 2005b). (A color version of this figure is available in the onUne journal.) 



pie. Note that the absolute normalizations of the AGN frac- 
tions in Fig. 15 are meaningless due to the above random- 
drawing scheme of ten galaxies vs. one AGN, and we are 
only interested in trends with color here. It seems clear, for 
a mass-matched sample, that there is no apparent clustering 
of AGNs in the CMDs and the AGN fraction does not vary 
much with color up to z ~ 2-3. We note that there appears 
to be a weak trend at 0.0 < z< 1.0 such that the fraction of 
X-ray luminous AGNs increases toward blue colors, which 
has also been seen in Sample C [see Fig. 10(D)]. This weak 
trend, which is unlikely to be due to AGN blue light con- 
tamination as demonstrated in § 4.6.3, implies a correspond- 
ing weak trend that the fraction of relatively low-luminosity 
AGNs (i.e., 41.9 < log(Lx/(erg s"')) < 43.0) increases toward 
red colors. These results confirm those obtained previously 
when considering Sample C. We repeated the construction 
of such a random mass-matched sample many times and find 
that these conclusions are stable. As a further check, we also 
constructed random non-mass-matched samples (i.e., for each 
AGN host, 10 unique galaxies with any stellar mass were ran- 
domly selected) and find that non-mass-matched samples do 
behave differently from mass-matched samples, i.e.: the re- 
sults obtained with non-mass-matched samples are different 
from those obtained with mass-matched samples, but are sim- 
ilar to those obtained with Sample B that is essentially non- 
mass-matched. This clearly demonstrates the mass-selection 
biases associated with non-mass-matched samples. 

Since the AGN fraction is approximately constant across 
colors for mass-matched samples, it is of interest to see 
whether the level of AGN activity is similar across colors. To 



this end, we plot in Fig. 16 histograms of X-ray luminosity (a 
proxy for the level of SMBH accretion) for AGNs with blue 
and red hosts in Sample B. We find that AGNs with blue and 
red hosts have similar distributions of X-ray luminosity (con- 
firmed by K-S tests), which indicates that AGNs with blue 
and red hosts have a comparable level of SMBH accretion (a 
similar result was obtained considering Sample C). 

The bottom panels of Fig. 15 show plots of SFR or upper 
limit on SFR as a function of redshift for AGN hosts and 
non-AGN galaxies in the above mass-matched sample. To 
make sample comparisons with upper limits included, we uti- 
lized the Astronomy SURVival Analysis package (ASURV; 
Isobe et al. 1986; Lavalley et al. 1992) to compute sam- 
ple means using the Kaplan-Meier (K-M) estimator and four 
nonparametric sample-comparison tests for censored data: the 
Gehan's generalized Wilcoxon test, the logrank test, the Peto 
& Peto generalized Wilcoxon test, and the Peto & Prentice 
generalized Wilcoxon test (these four tests typically produce 
similar results). We find, according to sample means derived 
using the K-M estimator as well as results from the above 
nonparametric sample-comparison tests, that (1) the SFRs of 
both AGN hosts and non-AGN galaxies generally increase to- 
ward higher redshifts from z ~ 0-3, which is consistent with 
the well-established formation epoch at z w 1-4; and (2) at 
< z < 1 the SFRs of AGN hosts are generally a factor of 
sa 2-3 [i.e., « 0.3-0.4 dex; see Fig. 15(c)] larger than those 
of non-AGN galaxies, while at 1 < z < 2 and 2 < z < 3 
AGN hosts and non-AGN galaxies on average have similar 
SFRs (i.e., the difference between their means of SFR is typ- 
ically within 0. 1-0.2 dex). The latter results are in agreement 
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Fig. 14. — (a - c) Color-mass diagrams for Sample B: AGN hosts are shown as filled circles and non-AGN galaxies are shown as small dots (colors have the 
same meaning as those in Fig. 10(a - i)). The two sets of dotted lines indicate the stellar-mass cut for Sample B (> IO'^Mq) and Sample C (> IO'^^Mq), 
respectively. Shown in the insets are the histograms of rest-frame U — V color for Sample B and Sample C, with AGN hosts shown as sohd lines and non-AGN 
galaxies as dotted lines (the peak value of each histogram has been rescaled to unity). (A - C) AGN fraction as a function of stellar mass in bins of Alog(Mt) = 0.5 
for Sample B (symbols have the same meaning as those in Fig. 10(A - 1)). (A color version of this figure is available in the online journal.) 



with previous results on AGN fraction vs. SFR obtained with 
Sample C in § 5.3. Our result that AGN hosts have elevated 
SFRs at < z < 1 also confirms the result from Silverman 
et al. (2009) who obtained a similar ratio between the SFRs 
of AGN hosts and non-AGN galaxies at < z < 1 [see their 
Fig. 9(f)]. 

We also find, if considering star-forming populations only 
[i.e., applying a cut of SFR > 10 (> 20, > 30) Mq yr"' 
to both AGN hosts and non-AGN galaxies at 0.0 < z < 1.0 
(1.0 < z < 2.0, 2.0 < z < 3.0)], that both AGN hosts and non- 
AGN galaxies have similar SFR distributions irregardless of 
redshift range considered. This result may contribute to ex- 
plaining the interesting trend we mentioned above that the 
factor of K, 2-3 difference between the SFRs of AGN hosts 
and non-AGN galaxies diminishes above z ~ 1 . The majority 
of AGN hosts/non-AGN galaxies at z > 1 have a SFR greater 

than w 5-10 M0 yr~', which means that their SFR distribu- 
tions are dominated by star-forming populations (thus being 
similar). Further investigation is required to understand the 
physical origin of AGN hosts having higher SFRs than non- 
AGN galaxies at low redshifts vs. AGN hosts and non-AGN 
galaxies having similar SFRs at high redshifts. 

There might appear to be an inconsistency between the two 
results obtained with mass-matched samples at < z < 1 that 
AGN hosts have higher SFRs than non-AGN galaxies and that 
the color distributions of AGN hosts and non-AGN galaxies 
are similar, since the former result suggests that AGN hosts 
should have bluer colors than non-AGN galaxies. However, 
as shown in § 5.3, the correlation between colors and SFRs 
is not very tight: for a given color, the SFRs span a wide 



range, although there may be a rough trend of bluer colors 
toward higher SFRs. Furthermore, the SFRs of AGN hosts 
are roughly a factor of 2-3 larger than the SFRs of non-AGN 
galaxies at < z < 1, which does not make a strong differ- 
ence in colors. Therefore, there is no real inconsistency in 
this issue. 

6. CONCLUSIONS AND SUMMARY 

In this work, we have assembled a sample of X-ray-selected 
moderate-luminosity AGNs as well as their parent galaxies 
in the most sensitive central areas of the 2 Ms CDFs to ex- 
tend color-magnitude relations of AGN hosts and non-AGN 
galaxies through the galaxy formation epoch. The uniqueness 
of this sample allows us to explore several important issues 
up to z ~ 2-3 when most of galaxy assembly occurred; such 
explorations at high redshifts have been largely unattainable 
in previous studies primarily due to data-depth issues. The 
main advantages of this work are (1) we have used the deep- 
est X-ray data, the best multiwavelength coverage, and spec- 
troscopic/photometric redshifts of the highest quality to date; 
and (2) we have used a template SED fitting method to derive 
galaxy physical properties (e.g., colors, stellar masses, and 
SFRs), which significantly reduces the failure rate of estima- 
tions for high-redshift sources whose SEDs are more prone to 
large uncertainties or errors. We summarize the main results 
of this work as follows: 

1. Non-AGN galaxy color bimodality (i.e., existence of 
the red sequence and the blue cloud) exists up to z ~ 3, 
with or without correcting for dust extinction. How- 
ever, there appears to be no apparent color bimodality 
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Fig. 15. — (Top, Middle) Same as Fig. 10, but for a mass-matched sample as defined in § 5.4. (Bottom) SFR or upper limit on SFR as a function of redshift 
for AGN hosts (shown as filled circles) and non-AGN galaxies (small dots) in the above mass-matched sample (symbols have the same meaning as those in 
Fig. 12(e)). The red (blue) crosses show the running means and associated errors of SFR (in bins of A^ = 0.5) for AGN hosts (non-AGN galaxies), which were 
computed using the Kaplan-Meier estimator provided by ASURV that fi'eats censored data (see § 5.4). The red crosses have been slightly shifted rightward for 
clarity of presentation. (A color version of this figure is available in the online journal.) 
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Fig. 16. — Histograms of X-ray luminosity for AGN hosts in Sample B. The red (blue) histograms are for red (blue) AGN hosts (red and blue AGN hosts are 
separated using Eqn. (4) in § 5.1.1). The two histograms in each panel are statistically indistinguishable. (A color version of this figure is available in the online 
journal.) 



for AGN hosts up to zw 2. These results (see § 5.1) 



confirm previous works and extend them to higher red- 
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shifts. 

2. For mass-matched samples up to z«2-3 (e.g., 
Sample C): (1) both AGN hosts and non-AGN galax- 
ies are luminous and share similar My distributions; (2) 
AGN clustering in the CMD is not apparent; and (3) 
the AGN fraction remains nearly constant (at ss 10%) 
regardless of color, which likely indicates the duty cy- 
cle of SMBH growth in typical massive galaxies (see 
§ 5.2 - § 5.4). These results are in contrast to those ob- 
tained with non-mass-matched samples up to z ss 2-3 
(e.g.. Sample B/D): (1) AGN hosts generally appear 
more luminous (i.e., having more negative values of 
My) than non-AGN galaxies; (2) AGNs predominantly 
reside in the red sequence, the top of the blue cloud, and 
the green valley in between (i.e., AGN clustering in the 
CMD); and (3) the AGN fraction generally increases as 
the color becomes redder (see § 5.2 and § 5.3). 

3. Luminous SMGs, the prototypical examples of extreme 
star-forming galaxies at z « 1-3, are spread throughout 
the CMD (see § 5.2). Likely due to dust-extinction ef- 
fects, they often do not show the blue colors typically 
associated with star-forming galaxies. This underscores 
the need for caution when studying the color-magnitude 
properties of small samples of remarkable sources. 

4. Most AGNs reside in massive hosts (i.e., M* > 
1O'"^M0), and the AGN fraction rises strongly toward 
higher stellar masses (see § 5.3), up to z ss 2-3, which 
confirms and extends previous works at lower redshifts. 

5. There are some strong trends of non-AGN galaxy /AGN 
host rest-frame colors, up to z~ 2-3 (see § 5.3): (1) 
the colors of both non-AGN galaxies and AGN hosts 
generally become bluer as the redshift increases; and 
(2) the colors of both non-AGN galaxies and AGN hosts 
generally become redder as the stellar mass increases. 

6. For mass-matched samples, the SFRs of AGN hosts are 
on average a factor of « 2-3 larger than those of non- 
AGN galaxies at z ~ 0-1, but this difference diminishes 

atz« 1-3 (see § 5.4). 

It has been demonstrated that galaxy evolution is best 
probed using stellar-mass selected samples (e.g., van Dokkum 
et al. 2006; Kriek et al. 2008; and references therein), 
van Dokkum et al. (2006) recommend that samples of high- 
redshift galaxies be selected by stellar mass, rather than by 
color or luminosity, because ( 1 ) the mass evolution of galaxies 
is probably gradual, whereas luminosities and colors can vary 
dramatically on short timescales due to starbursts and dust, 
and (2) models of galaxy formation can predict masses with 
higher confidence than luminosities and colors. Moreover, 
stellar-mass selected samples make comparisons between the 
properties of AGN hosts and non-AGN galaxies more appro- 
priate. For example, Silverman et al. (2009) showed that 
a mass-selected sample is required to reduce the strength of 
an artificial peak in the AGN fraction falling in the transi- 
tion region (i.e., the green valley) due to the fact that many 
blue-cloud galaxies have low mass-to-light ratios in luminos- 
ity limited samples. 

In this work, we have further demonstrated the importance 
of using mass-matched samples for investigations of color- 
magnitude relations of active and non-active galaxies, because 



results and conclusions obtained with a mass-matched sample 
are different from those obtained with a non-mass-matched 
sample (see § 5.2 - § 5.4). For instance, for a non-mass- 
matched sample, AGNs predominantly reside in the red se- 
quence, the top of the blue cloud, and the green valley in 
between; this has been presented as evidence for AGN feed- 
back because the location of AGN hosts in the CMD is con- 
sistent with the transitional region where the blue-to-red mi- 
gration of galaxies occurs due to the quenching of star for- 
mation (e.g., Nandra et al. 2007). However, the AGN clus- 
tering in the CMD largely disappears and the AGN fraction 
remains nearly constant as a function of color when consid- 
ering a mass-matched sample. Therefore, in order to obtain a 
complete picture about the role of AGNs in galaxy evolution 
by means of CMD analyses, one must fully take into account 
the dependence of AGN activity on stellar mass (i.e., using 
mass-matched samples). 

With mass-selection effects taken into account, we find that 
the main results in this work (except for color-bimodality, 
which likely needs some quenching mechanism) can be rea- 
sonably explained by a combination of two main ingredients: 
(1) AGNs at z~ 0-3 preferentially reside in massive galax- 
ies that generally tend to have redder colors (i.e., color-mass 
correlation); and (2) galaxies evolve passively (e.g., Kodama 
& Arimoto 1997; Stanford, Eisenhardt, & Dickinson 1998; 
van Dokkum & Franx 2001; Daddi et al. 2004; Bothwell, 
Kennicutt, & Lee 2009) or secularly (e.g., Kormendy & Ken- 
nicutt 2004; Hopkins & Hernquist 2006); this ingredient is 
mainly responsible for the evolution of galaxies, i.e., blue- 
to-red migration. Consequently, our results tightly constrain 
any effects from moderate-luminosity AGN feedback upon 
the color-magnitude properties of galaxies fromz ~ 0-3 when 
most of galaxy assembly occurred. We note that moderate- 
luminosity AGNs, like those studied here, dominate the ac- 
cretion density and SMBH growth only at z ~ 0-1 while more 
luminous AGNs dominate above z ~ 1 (e.g., Ueda et al. 2003; 
Hasinger, Miyaji, & Schmidt 2005). For moderate-luminosity 
AGNs, some simulations predict that feedback may not be 
important (e.g., Hopkins & Hernquist 2006), whereas other 
models suggest it is (e.g., Fabian, Vasudevan, & Gandhi 2008; 
Raimundo et al. 2010). Given that such simulations and mod- 
els generally rely on significant assumptions and are there- 
fore subject to large uncertainties, further observational tests 
of moderate-luminosity AGN feedback are required. 

We show that the AGN fraction is approximately constant 
across colors for galaxies with similar stellar masses and that 
AGNs with blue and red hosts have comparable X-ray lumi- 
nosities (see § 5.4). These results imply that AGN activity 
is about as prevalent in massive blue galaxies as it is in red 
galaxies. Assuming that the SMBH mass generally decreases 
toward blue colors since the bulges may be smaller in blue 
(late-type) hosts than in red (early-type) hosts for a given 
galaxy stellar mass, the less-massive SMBHs in blue hosts 
may undergo more significant fractional growth (i.e., with 
larger Eddington ratios) than those in red hosts. If this effect 
exists and is significant, AGN feedback may be more likely to 
occur in those systems having higher ratios of LAGNZ-'Wbuige- 
However, the above assumption is not certain. For example, 
Schawinski et al. (2010) show that in local universe there 
is no simple relationship between colors and morphological 
types [e.g., a red (blue) galaxy is not necessarily an early-type 
(late-type) galaxy and vice versa]; they also show that both 
early- and late-type AGN hosts have similar typical SMBH 
masses and Eddington ratios (although both SMBH masses 
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and Eddington ratios span a broad range). At high redshifts, it 
is more challenging to obtain reliable morphological informa- 
tion and estimates of SMBH masses. Therefore, it is difficult 
to assess how significant AGN feedback effects could be in 
some massive blue galaxies if such feedback even exists at 
z ~ 0-3, given the current constraints. Further investigation 
of these effects would be worthwhile. 

A great deal of further work, in addition to that already 
mentioned, can be done to improve and extend the analyses 
presented here, e.g., obtaining even deeper X-ray and near- 
infrared data. According to Fig. 3(a), we would expect to 
identify an AGN with Lx > 10"*^ erg s"' out to z « 3, but we 
would not expect to identify an AGN with Lx a few times 
10^^ erg s"' beyond z « 1.5. We note that this is already the 
best that the currently deepest X-ray surveys, the 2 Ms CDFs, 
can offer. Therefore, X-ray data with even higher sensitiv- 
ity are required to provide better completeness, which will 
also increase the size of the AGN sample by detecting more 
highly obscured and moderate-luminosity AGNs at high red- 
shifts. Obtaining ultradeep near-infrared data would also be 
beneficial for the analyses here. First, ultradeep near-infrared 
selected source catalogs would include a population of less- 
massive red galaxies and thus lower the stellar-mass com- 
pleteness limit. Second, ultradeep near-infrared data have 
proved critical to improve Zzphot estimates and thus estimates 



of other source properties, especially for high-z sources, e.g., 
Wardlow et al. (2010) made use of ultradeep HAWK-I J 
and K^ data as well as other multiwavelength data to derive 
high-quality photometric redshifts for a sample of E-CDF-S 
SMGs. Another valuable project within the framework of 
color-magnitude and color-mass relations, as an extension of 
our work here, would be to constrain the colors of quasar hosts 
out to high redshifts, although practically this is difficult be- 
cause the host light is overwhelmed by quasar emission except 
for highly obscured quasars. 
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